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FOREWORD 


TMs  is  t.hp  final  report  covering  Phase  II,  Parametric  and  Systems 
Studies,  of  the  work  performed  under  the  Department  of  the  Navy,  Bureau  of 
Naval  Weapons,  Contract  NOw  62-0887-d,  in  accordance  with  the  Contractor's 
proposal  "Jet  Noise  Suppression  Research  and  System  Studies",  P-61-86,  dated 
March  1961. 

The  objectives  of  Phase  II  of  the  program  were  t©  evaluate  the  effects  of 
noise  from  aircraft  carrier  deck  Jet  engine  operation  oa  hearing,  speech  inter 
ference,  and  structures,  for  the  purpose  of  establishing  suitable  acoustical 
criteria  for  in-flight  jet  noise  suppressors .  In  addition,  the  effects  of 
suppressor  aerodynamic  performance  and  weight  on  aircraft  mission  performance 
were  evaluated  parametrically  for  several  aircraft  missions. 

This  report  is  presented  in  two  volumes .  Volume  I  is  concerned  with  the 
effects  of  jet  engine  noise  on  hearing,  speech  interference,  and  sonic  fatigue 
and  Volume  II  deals  with  the  effects  of  noise  suppressors  on  aircraft  mission 
performance . 


I .  INTRODUCTION 


This  report  is  part  of  a  study  undertaken  for  the  Navy  in  an  attempt  to 
better  define  the  parameters  relative  to  the  problem  of  an  in-flight  noise 
suppressor  for  carrier  based  jet  aircraft. 

The  design  of  such  a  suppressor  must  represent  the  result  of  compromise. 
The  object  is  to  balance  acoustical  performance  of  the  suppressor  against 
factors  which  influence  aircraft  performance  -  such  as  internal  pressure  loss, 
external  drag  loss ,  and  weight  -  and  also  against  any  ill  effects  which  noise 
may  have  on  the  aircraft  or  on  carrier  personnel. 

Part  of  the  overall  problem,  then,  is  to  study  all  possible  effects  of 
the  high  noise  levels  characteristic  of  aircraft  operations,  and,  where 
possible,  to  recommend  criteria  relative  to  noise  levels,  exposure  time,  etc. 
This  divides  naturally  into  two  categories;  first,  the  effect  of  noise  on 
people,  and  second,  the  effect  of  noise  on  structure. 

In  the  first  category,  the  two  major  problems  are  hearing  damage  risk, 
and  interference  with  communication.  Although  much  has  been  written  on  both 
of  these  fields,  there  is  need  to  coordinate  the  available  information,  and 
to  present  it  in  a  form  more  directly  applicable  to  the  problem  at  hand. 

A  similar  statement  might  also  be  made  relative  to  the  effect  of  noise  on 
structures . 

This  report  is  based  almost  entirely  on  literature  reviews,  and  discus¬ 
sions  held  with  research  workers  in  the  several  fields.  No  tests  were  con¬ 
ducted  as  part  of  tne  project. 
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II.  SUMMARY  AND  CONCLUSIONS 


Inf o rma t ion  gathered  from  a  literature  search  -  notably  the  work  of 
3*  7 

Ward,  Ctorig,  et  al  -  relative  to  hearing  loss  vs.  noise  exposure  is 
reviewed.  The  results  of  such  studies  are  applied  to  the  specific  problem 
of  determining  criteria  applicable  to  the  hearing  damage  risk  problem  aboard 
aircraft  carriers.  Simple  nomograms  {Figs.  4,  8,  and  9)  are  derived  to 
quickly  evaluate  the  maximum  allowable  exposure  time  to  a  given  sequence 
of  take-offs  expressed  in  terms  of  octave-band  sound  level  vs,  time  plots 
at  the  observer's  ear.  These  charts  are  based  upon  the  allowance  of  a  com¬ 
puted  TTS^  (temporary  threshold  3hift  2  minutes  after  exposure  to  the  noise) 
of  10-13  db  at  2000  cps  and  20-23  db  at  4000  cps.  Additional  nomograms  allow 
rapid  computation  of  TTS  growth  and  recovery  during  noisy  and  quiet  periods 
respective  ly . 


TTS  is  considered  to  be  of  prime  importance  in  specifying  noise  criteria 

since  studies^  indicate  that  a  noise  exposure  which  results  in  a  daily 

temporary  threshold  shift  of  X  db  should,  in  a  10  year  period,  produce  a 

permanent  threshold  shift  of  X  db. 

14  9 

The  Noise  Cumulator  *  is  discussed  briefly,  and  an  example  is  given 
of  the  application  of  such  a  device  to  TTS  predation. 

Speech  communication  In  the  presence  of  noise  is  discussed  briefly  in 
terms  of  the  well-known  speech- inter ference- level  (SIL) ^  and  also  in 

20  2  1 

terms  of  the  more  fundamental  20  band  Al  (Articulation  Index)  calculation.  J 

27 

A  simplified  6  octave-band  Al  calculation  recently  proposed  by  Xryter  ,  is 
discussed  in  some  detail.  This  approach  has  the  advantage  of  greater  accuracy 
than  can  be  obtained  by  the  use  of  S1L  with  but  a  minimum  of  added  complica¬ 
tion.  Furthermore,  modi f i cat  1 ons  to  the  Al  calculation  may  be  employed  to 
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take  into  account  various  influencing  factors  which  are  completely  neglected 
with  the  usual  SIL  approach. 

Several  examples  are  presented  which  illustrate  the  simplified  6  band 
AI  calculation,  and  also  serve  to  indicate  its  superiority  over  normal  SIL 
usage . 

The  effect  of  noise  on  structures  is  examined  in  several  phases.  In 
general,  we  follow  the  approach  of  Fitch^  which  considers  first  the  determin¬ 
ation  of  the  acoustic  environment,  and  the  duration  of  important  noise  expos¬ 
ures  as  based  on  a  typical  mission  analysis.  The  subject  of  structural  design 
criteria  is  then  approached  via  accelerated  discrete  frequency  life  testing, 
the  results  of  which  are  interpreted  using  a  sine-random  equivalence  analysis. 
The  theory  of  vibration  excitation  due  to  random  noise  and  the  mechanism  of 
structural  fatigue  are  discussed. 

Design  charts  developed  by  Belcher4^  and  McGowen^  are  discussed.  These 
charts  relate  sound  pressure,  fatigue  life,  resonant  frequency,  and  dimensions 
of  several  common  panel  constructions. 

Because  of  the  very  nature  of  the  problem,  it  is  not  possible  to  specify 
sound  criteria  as  fixed  numbers  applicable  to  all  situations.  However,  the 
principles  and  information  presented  in  this  report  will  allow  the  calculation 
of  maximum  recommended  exposure  time  to  specific  aircraft  and  flight  operation 
with  personnel  at  specified  locations.  They  will  also  serve  to  specify  design 
parameters  for  speech  communication  systems  for  given  situations.  In  addition, 
they  allow  a  rapid  prediction  of  the  increase  in  allowable  personnel  exposure 
time,  improvement  of  speech  communication,  and  reduction  in  danger  of  structural 
fatigue  failures  due  to  the  use  of  a  given  exhaust  silencer. 

It  is  proposed  therefore,  that  detailed  studies  be  undertaken  to 
evaluate  the  above  aspects  for  all  situations  of  major  interest. 
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III.  EFFECT  OF  NOISE  ON  PEOPLE 


1.  Hearing  Damage  Risk  Criteria 

Years  ago,  noise  criteria  relative  to  hearing  damage  risk  were 

expressed  simply  as  weighted  sound  levels,  octave-band  levels,  sone 

values,  etc.  that  should  not  be  exceeded.  It  became  apparent,  however, 

that  exposure  time  is  an  important  factor  which  cannot  be  neglected. 

Several  methods  of  setting  criteria  involving  both  levels  and  exposure 

time  were  then  proposed  -  the  one  most  widely  adopted  being  the  equal 

energy  hypothesis  often  expressed  as 

IT  =  constant  (1) 

where  I  represents  sound  intensity  and  T  represents  time.  ^  The  Air  Force 

2 

Regulation  160-3  (1956)  is  based  upon  this  assumption. 

The  above  equation  states  that  if  the  intensity  of  noise  (perhaps 
in  a  given  frequency- band)  is  doubled  (increased  by  3  db) ,  then  only 
half  the  previous  exposure  time  should  be  permitted.  Similarly,  if  the 
noise  is  increased  in  intensity  by  a  factor  of  100  (20  db) ,  then  the 
exposure  time  should  be  reduced  to  17*  of  the  original  time.  With  the 
high  noise  levels  typical  of  modern  jet  engines,  such  a  process  gives 
rise  to  extremely  short  allowable  exposure  times  -  or  alternatively  - 
more  required  ear  protection  than  is  available  from  standard  ear  plugs 
and  muffs. 

In  recent  years,  both  studies  of  controlled  noise  exposure  vs. 

3-  3 

hearing  loss  and  studies  of  hearing  losses  actually  suffered  by 

i  4  .  .  9 , 10  .  , . 

people  exposed  to  jet  engine  noise  indicate  that  the  equal  energy 

hypothesis  is  somewhat  overprotec tive  for  high  level  noise.  This  is 
particularly  true  for  intermittent  bursts  of  intense  noise^  -  for  which 
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Eq.  (1)  may  be  rewritten  as 

IRT  «  constant  (2) 

where  R  is  defined  as  the  "on- fraction'*  of  the  noise. 

g 

One  of  the  most  recent  papers  in  this  field  is  that  of  Kryter. 
However,  although  he  gives  damage  risk  criteria  for  exposure  to  continu¬ 
ous  noise,  his  work  is  not  directly  applicable  to  the  problem  of  inter¬ 
mittent  noise  -  which  is,  of  course,  of  concern  during  launch  operations 

aboard  an  aircraft  carrier.  Fortunately,  though,  the  work  of  Ward, 

3-  7 

Glorig,  et  al  is  applicable  to  this  problem  and  may  be  used  to  estimate 
both  temporary  and  permanent  hearing  threshold  shifts  when  one  is  exposed 
to  a  complex  time-varying  noise  such  as  that  produced  by  a  sequence  of 
jet  aircraft  take-offs  from  a  carrier  deck. 

Before  considering  this  specific  problem,  however,  we  will  first 
summarize  the  basic  findings  of  these  studies.  First,  the  general 
equation  relating  (temporary  threshold  shift  two  minutes  after 

cessation  of  noise)  to  a  noise  exposure  of  duration  T  minutes  at  a  sound 
pressure  level  S  is  apparently  of  the  form: 

TTS2  =  K(S  -  So)  [log10  (T/To)]  +  C  (3) 

4 

where  K,  T^,  and  C  are  constants  for  a  given  type  of  noise.  (Two 

minutes  is  usually  taken  as  a  suitable  time  tor  measurement  since  hearing 

recovery  is  erratic  during  the  first  minute  or  so  after  the  noise  ceases.) 

Second,  TTS  at  a  given  frequency  is  primarily  determined  by  noise  in 

the  octave- band  near  that  frequency  and  in  the  octave- band  immediately 

below  this  one.  In  other  words,  4  KC  loss  is  largely  due  to  noise  in  the 

4 

1200-2400  and  2400-4800  cps  bands,  etc.  The  following  tabulation  gives 
empirical  equations  for  growth  of  TTS^  hue  to  octave-band  noise.  In  these 

equations  it  is  assumed  that  T  is  of  the  order  of  5  minutes  or  more. 
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TABLE  I 


Exposure 

Test 

Band 

Frequency 

Equation  for  TTS2 

600-1200 

1500 

0. 53(S-71)  (log  T  -  0.44)  -  3 

600-1200 

2000 

0. 41(S-68)  (log  T  +  0.  15)  -  8 

1200-2400 

3000 

0. 38(S-65)  (log  T  +  0.55)  -  13.5 

1200-2400 

4000 

0.61(S-70)  (log  T  +  0.33)  -  9.5 

2400-4800 

4000 

0. 91(S- 75)  (log  T  +  0.  19)  -  8 

2400-4800 

6000 

0. 51(S-68)  (log  T  +  1.80)  -  22 

Although  the  tests  were  confined  to  octave- band  levels  of  120  db  and  below, 
it  is  not  unreasonable  to  extrapolate  the  above  equations  to  octave-band 
levels  approaching  130  db.  However,  above  130  db,  even  very  short 
exposures  are  considered  serious.  (Also,  exposures  for  which  the  expected 
TTS2  exceeds  30  db  should  not  be  permitted), 

Third,  if  the  noise  is  intermittent  with  bursts  of  duration  between 
1/4  second  and  1  minute,  the  computed  TTS2  as  given  by  equation  (3)  should 
be  multiplied  by  R,  the  "on- fraction"  of  the  noise, ^  Moreover,  if  the 
noise  consists  of  several  discrete  levels  with  different  on- fractions , 
it  appears  that  the  equation  should  be  written; 

n 

TTS.  =  K  [  R.  (S  ,-S  )  log(T/T  )  +  C  (4) 

Z  1X0  o 

Fourth,  after  cessation  of  noise,  recovery  at  all  frequencies  may 
be  approximated  by  the  relation 

TTSC  =  (TTS2  +9)  (1-0.27  log10  i)  -  9  (5) 

4 

where  the  recovery  time  t  is  assumed  to  be  in  excess  of  2  minutes. 

Fifth,  if  repeated  exposures  to  noise  are  separated  by  periods  of 
relative  quiet  with  durations  exceeding  about  1  minute,  it  is  necessary 
to  employ  a  cumulative  calculation  for  TTS  growth  and  recovery. ^  The 
specific  steps  are  outlined  below. 
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1. 


Compute  TTS^  after  first  exposure* 

2.  Compute  TTS^  after  t  minutes  recovery. 

3.  Using  the  growth  equation,  calculate  the  value  of  T  that  would 
be  required  to  obtain  the  TTS^.  of  step  (2)  for  a  TTS^* 

4.  Add  this  "equivalent  exposure  time"  to  the  actual  T  for  the 
second  exposure  and  compute  TTS2  for  the  total. 

3.  Repeat  this  process  for  the  entire  sequence  of  exposures  and 
recoveries . 

Sixth,  if  more  than  one  octave-band  contributes  appreciably  to  the 
computed  TTS^,  an  approximate  combination  of  effects  may  be  made  by  adding 
the  expected  TTS  values  as  one  normally  adds  decibel  quantities.  That  is, 
with  a  TTS  2  of  10  db  at  2000  cps  due  to  600-1200  cps  noise  and  another 
TTS2  of  10  db  at  2000  cps  due  to  1200-2400  cps  noise,  one  would  expect  a 
total  TTS  of  13  db  at  2000  cps. 

Seventh,  we  have  the  question  of  permanent  threshold  shift  (PTS)  versus 
temporary  threshold  shift  (TTS).  Studies  of  various  industrial  environ¬ 
ments  indicate  that  dally  exposures  producing  X  db  TTS  will  after  10-13 
years  give  rise  to  a  PTS  of  the  order  of  X  db.^  At  4  KC,  this  PTS  is 
approximately  a  maximum.  That  is,  further  exposure  does  not  increase  it 
appreciably.  For  2  KC  the  PTS  continues  to  increase  with  further  exposure. 

A  somewhat  simplified  analysis  of  the  above  type  is  included  in  the 
Proposed  International  Standard  ISO  43  (Secretariat- 194) 314E ,  June  1963 
for  ,rNoise  Rating  with  Respect  to  Conservation  of  Hearing,  Speech  Communi¬ 
cation,  and  Annoyance".  The  curves  and  tabulations  in  this  document  are 
based  largely  on  an  allowable  TTS 2  of  12  db  at  2000  cps  for  a  daily 

exposure.  Kryter  suggests  10  db  at  and  below  1000  cps,  13  db  at  2000  cps, 

8 

and  20  db  at  and  above  3000  cps.  With  losses  appreciably  greater  than 
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these  figures,  ability  to  understand  conversation  becomes  severely  impaired* 
The  proposed  ISO  method  is  not  directly  applicable  to  our  problem  as 
it  does  not  include  the  R  factor  approach  for  short  time  variations  in 
exposure.  Moreover,  a  simplified  computation  is  made  only  for  the  one 
frequency  band  which  is  highest  with  respect  to  a  specific  set  of  rating 
curves.  The  procedure  explicitly  outlined  above  should  be  somewhat  more 
accurate  in  the  latter  regard. 

In  the  interests  of  simplicity,  however,  it  seems  sufficient  to 
consider  TTS  only  at  2000  and  4000  cps .  Actually,  PTS  at  4000  cps  is  not 
very  important  from  the  standpoint  of  understanding  conversation,  and  at 
1000  cps  the  TTS  will  almost  invariably  be  considerably  less  than  at  2000 
cps.  In  a  recent  private  communication.  Dr.  Glorig  3tated  that  the  second 
equation  of  Table  I  might  also  he  used  ror  TTS^  at  2000  cps  due  to  noise 
in  the  -200-2400  cps  band.  Let  us  consider,  then,  the  following  four 
equations : 

TABLE  II 


Exposure 

Band 


Test 

Frequency 


Equation  for  TTS^ 


600-1200 

2000 

1200-2400 

2000 

1200-2400 

4000 

2400-4800 

4000 

0. 4l(S- 68)  ( log  T  +  0. 15)  -  8 

0. 41(S- 68)  (log  T  +  0. 15)  -  8 

0. 61(S-  70)  (log  T  +  0.33)  -  9.5 

0. 9 1(S- 75)  (log  T  +  0.  19)  -  8 


These  equations  assume  the  noise  to  be  broadband  in  character.  With  pure 

tones  (such  as  from  compressor  noise)  the  problem  is  more  serious.  The 

8 

work  of  Kryter  indicates  that  in  such  cases  the  above  equations  should 
be  modified  by  subtracting  5  db  from  the  various  values  given  for  - 
since  the  critical  bandwidth  for  TTS  appears  to  be  of  the  order  of  1/3 
octave.  Alternatively,  the  equations  may  be  used  as  given,  and  5  db 
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added  to  the  values  employed  for  S. 

We  will  now  attempt  to  illustrate  the  above  approach  by  computing 
the  expected  TTS^  at  2000  and  4000  cps  for  a  man  wearing  ear  plugs  while 
exposed  to  aircraft  take-off  operations  on  a  carrier  deck.  (Ref.  11 
indicates  that  the  noise  from  landing  operations  is  less  important). 

Fig.  1  indicates  a  typical  idealized  exposure  time  pattern  with  after- 

u  ,  11,12 

burner  operation. 


time- seconds 
Fig.  1 


Here  alternate  take-offs  would  probably  be  made  from  two  catapults,  so 
for  simplicity  we  will  assume  that  the  man  is  halfway  between  the  two, 
at  a  distance  of  25  feet  from  the  tail  of  each  aircraft,  and  at  the  angle 
of  maximum  noise  radiation.  The  level  represents  military  operation, 
wnile  S2  indicates  the  increased  levels  when  afterburners  are  started. 
Compressor  noise  is  considered  to  be  negligible  with  respect  to  exhaust 
noise,  so  the  equations  of  Table  II  will  be  used  without  any  5  db  correc¬ 
tions.  It  will  be  assumed  that  20  aircraft  take  off,  giving  a  total 
exposure  time: 

T  =  0.25  +  19  x  0.50  -  9.75  minutes 

Four  such  flights  will  be  repeated  at  2  hour  intervals,  and  TTS  computed 
for  the  entire  operation.  It  is  also  assumed  that  all  aircraft  are  North 
American  A3J*lfs,  each  using  two  General  Electric  J79  engines. 
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From  Ref.  13,  the  sound  levels  to  be  expected  are  given  in  Table  III  for 
both  military  and  A/B  operation.  Also  tabulated  is  the  "average"  attenua¬ 
tion  for  V-51R  ear  plugs,  as  taken  from  Ref.  10.  Subtracting  these  values 
gives  the  estimated  exposure  for  the  two  periods.  Next,  the  effective 
"on- factor”  is  computed,  and  finally  the  TTS^  after  the  first  group  of 
planes  take  off.  These  are  found  to  be  7.2  and  12  db  at  2000  and  4000  cps, 
respectively  -  after  combining  the  individual  components  as  previously 
outlined.  At  the  bottom  of  Table  III,  an  average  growth  equation  is  given 
for  this  combined  TTS^. 

Next,  we  consider  recovery  during  the  following  110  minutes  of  assumed 
quiet.  (Other  noise  could,  of  course,  effect  this  recovery).  For  t  =  110, 
equation  (5)  becomes 


TTS110  «  0. 53  TTS2  -  4.2 


As  shown  in  Table  IV,  recovery  is  essentially  complete  at  2000  cps,  but 
there  is  a  residual  TTS  of  2.1  db  at  4000  cps.  This  2.1  db  is  then  used 
with  the  average  growth  equation  of  Table  III,  to  compute  an  equivalent 
exposure  of  2.4  minutes.  Hence,  the  effective  T  for  the  second  group  of 
planes  Is  9.8  +  2,4  =  12.2  minutes  -  which  from  the  growth  equation  should 
produce  13.6  db  TTS^  at  4000  cps.  This  process  may  then  be  repeated  for 
the  third  and  fourth  flights  as  indicated. 


TABLE  III 


Octave-Band 

si 

S2 

V-51R  Attenuation 

(with  plugs) 

S2  (with  plugs) 


2000  cps 


600 

1200 

1200 

2400 

138 

140 

144 

147 

-18 

-25 

120 

115 

126 

117 

4000  cps 


1200 

2400 

2400 

4800 

140 

133 

147 

140 

-25 

-30 

115 

103 

117 

110 
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TABLE  III  (Cont'd) 


2000  cps 

4000  cps 

r  <s  -s  ) 

Rj(S^-5°) 

1/4  x  52 

1/4  x  47 

1/4 

x  45 

1/4  x  28 

1/4  x  58 

1/4  x  49 

1/4 

x  47 

1/4  x  35 

LRi(si"V 

27.5 

24 

23 

15.75 

(log  X  +  log  To) 

1.14 

1.14 

1.32 

1.18 

K 

0.41 

0.41 

0.61 

0.91 

TTS2 

4.9 

3.3 

9.0 

8.9 

Combined  TTS2 

7.2 

12 

Aver*  Growth  Eqn. 

13. 3( log  T 

+  0.15)  -  8 

16. 4( log  T 

+  0.26)  -  8.5 

TABLE  IV 

2000  cps 

4000  cps 

1st  TTS2 

7.2 

12 

118  no 

— 

2.1 

Equiv.  T. 

2*4 

Total  T 

12.2 

2nd  TTS2 

7.2 

13.6 

TTS110 

— 

3.0 

Equlv.  T. 

2.8 

Total  T 

12.6 

3rd  TTS2 

7.2 

13.8 

TO110 

— 

3.1 

Equiv.  T 

2,8 

Total  T 

12.6 

4th  TTS. 

7.2 

13.8 
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Rather  than  calculate  as  indicated  in  Table  III,  one  might 

compute  the  numerical  average  of  the  sound  levels  3^  and  S2  (since  in  this 
example  they  persist  for  equal  time  intervals)  and  use  an  R  factor  of  0*5 
(since  the  periods  of  noise  and  relative  quiet  are  equal).  Thus,  in  the 
first  column  of  Table  III,  we  might  compute  Sav  *»  123,  and 


TTS2  "  0.41  x  0.5  (123-68)  1.14  -  8 
*  4.9 

Actually,  this  procedure  is  identical  in  principal  to  the  one  employed 
above;  but  is  easier  to  use.  It  may  be  readily  extended  to  more  complex 
sound  level  histories  by  calculating  the  weighted  numerical  average  of 
the  levels  -  possibly  with  the  aid  of  a  planimeter.  (Note  that  we  are 
not  computing  the  level  of  the  mean  sound  pressure,  or  mean- squared 
pressure).  Thus,  for  the  case  illustrated  in  Fig,  2 


Fig.  2 


the  average  level  is 

Sav  -  1/7(115  +  115  +  120  +  130  +  125  +  125  +  125) 
*  122  db 
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as  Indicated  by  the  dotted  line,  and  a  factor 

v1 

R  =  ~  0.637 

is  applicable  to  the  term  (Sav  -  3q) 

When  the  growth  equations  are  employed  in  this  manner,  it  becomes 
possible  to  put  them  in  the  form  of  simple  nomograms.  Fig.  3  is  such 
a  nomogram  applicable  to  both  the  600-1200  cps  band  and  the  1200-2400  cps 
band  tor  estimating  TTS^  at  2000  cps.  Thus,  for  the  calculation  in  the 
first  column  of  Table  III,  we  draw  a  line  between  S&v  *  123  and  R  «  0.5, 
note  the  intersection  on  the  unlabeled  turn  line,  and  draw  a  line  from 

this  point  to  T  =  9.8  intersecting  the  TTS^  scale  at  5  db. 

Now,  if  as  a  criterion,  we  do  not  allow  the  TTS2  at  2000  cps  as 
computed  from  Fig.  3  to  exceed  10  db  for  either  the  600-1200  band  or  the 

1200-2400  band,  then  the  total  TTS^  at  2000  cps  should  not  exceed  13  db. 

This  figure  is  consistent  with  present  day  thinking  -  as  previously  noted. 
With  this  10  db  criterion.  Fig.  3  becomes  very  useful  for  estimating 
allowable  combinations  of  S&v,  R,  and  T.  For  example,  assuming  R  *  0.5 
we  may  tabulate  as  follows: 

TABLE  V 

S  Max.  allowable  T  (minutes) 

flv 

125  25 

120  35 

115  52 

For  TTS2  c  10,  Fig.  3  may  be  simplified  to  Fig.  4  which  permits  even  quicker 
computations  of  the  above  type  -  though  with  unavoidable  loss  in  perspec¬ 
tive,  with  regard  to  the  relative  variation  of  the  factors  involved. 

(For  R  «  1,  Fig.  3  reduces  to  a  criterion  for  continuous  noise;  and  though 
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still  applicable  it  could  obviously  be  simplified  to  a  3  line  chart)* 

Where  repeated  exposures  at  intervals  are  of  concern  (as  in  the 
flight  deck  example  previously  discussed)  we  must  consider  TTS  recovery 
between  exposures.  Fig.  5  was  constructed  from  Eq.  (3)  to  make  such  a 
computation  easy.  As  a  suitable  criterion  here,  it  is  suggested  that  the 
time  interval  between  such  exposures  be  chosen  with  the  aid  of  Fig,  5  so 
that  the  computed  TTSt  does  not  exceed  1  db.  In  that  event,  exposures 
which  produce  the  maximum  TTS2  of  10  db  should  not  be  repeated  more 
frequently  than  once  every  two  hours,  though  exposures  to  less  intense 
noise  could  be  more  frequent. 

If  desired,  Figs.  3  and  5  may  be  employed  to  determine  the  expected 
TTS  after  a  ser'es  of  repeated  exposures  -  as  in  the  tabulated  example. 

If  only  one  of  the  two  octave-bands  is  of  importance,  the  method  of 
computation  appears  obvious.  However,  if  neither  band  is  negligible,  the 
use  of  an  average  growth  equation  may  be  effected  by  using  Fig.  3  in  a 
s lightly  dif ferent; manner .  Suppose,  for  example,  that  we  found  TTS2  3  9  db 
for  both  the  600-1200  band  and  the  1200-2400  band.  We  would  then  assume 
the  total  TTS2  to  be  12  db  at  2000  cps.  Now,  taking  the  exposure  time, 

T,  as  10  minutes,  we  may  draw  a  line  through  T  =  10  and  TTS2  89  12  and 
determine  the  intersection  on  the  unmarked  turn  line,  (This  is  the 
effective  combination  of  R  and  S  ,  but  need  not  be  evaluated  numerically.) 
Suppose  next  that  we  desire  a  1  hour  interval  between  repeated  exposures. 
How  does  TTS  build  up  over  several  9uch  cycles?  From  Fig.  5,  we  find 
TTSt  -  4  db  at  1  hour.  Inserting  this  value  on  the  TTS2  scale  of  Fig.  3, 
and  using  the  turn  line  point  previously  described,  we  find  an  equivalent 
exposure  time  of  about  3.4  minutes.  Hence,  the  effective  time  for  the 
second  exposure  is  13.4  minutes.  Using  this  same  turn  line  point,  we 
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now  finds  TTS2  =  14  db.  Similarly,  after  another  cycle  TTSj  *  15  db,  etc. 

Figs.  6  and  7  are  nomograms  giving  the  same  type  of  TTS2  calculation 

for  lost  at  4000  cps  due  to  1200-2400  cps  noise  and  2400- 4b00  cps  noise, 

g 

respectively.  Here,  in  line  with  Kryter's  reasoning  ,  we  suggest  an 
allowable  maximum  of  20  db  as  computed  for  either  band.  If  only  jet 
exhaust  noise  were  of  concern,  it  is  probable  that  the  criterion  of  Fig.  3 
would  govern  the  problem.  However,  compressor  whine  may  also  be  of 
importance,  and  may  easily  be  sufficient  to  require  the  use  of  Figs.  6  and 
7.  (As  noted  previously,  5  db  should  be  added  to  Sflv  for  the  case  of 
pure  tones.)  Since  recovery  seems  to  be  independent  of  frequency,  Fig.  5 
may  still  be  used  at  4000  cps.  Also,  as  before,  it  is  possible  to  effec¬ 
tively  use  an  average  growth  equation  if  both  octave  bands  are  of  importance 
For  such  a  calculation,  one  may  use  either  Fig.  6  or  Fig.  7  with  Fig.  5  - 
the  difference  in  the  final  result  being  negligible. 

Figs.  8  and  9  correspond  to  Fig.  4,  except  that  these  are  drawn  for 
a  limiting  TTS2  of  20  db  at  4000  cps  due  tc  noise  in  the  1200-2400  and 
2400-4800  bands. 

In  summary,  to  use  this  approach  in  its  simplest  form  as  a  criterion 
to  determine  the  maximum  allowable  exposure  time  to  a  series  of  take-offs 
as  described  previously,  one  need  only  determine  plots  similar  to  Fig.  2 
(evaluating  both  S  and  R)  for  each  of  the  3  bands  600-1200,  1200-2400, 
and  2400-4800  cps.  Then  Fig.  4  should  be  used  with  both  the  600-1200  and 
1200-2400  cps  bands.  Fig.  8  with  the  1200-2400  band  only,  and  Fig.  9  with 
the  2400-4800  band,  with  5  db  additions  to  S  for  any  strong  pure  tones. 

The  minimum  value  for  T  obtained  among  these  4  computations  represents  the 
maximum  allowable  exposure.  If  desired,  of  course,  more  detailed  informa- 
tion  on  expected  TTS  growth  and  recovery  can  be  obtained  from  the  other 


charts  as  described  above. 
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In  most  cases,  personnel  will  be  wearing  some  form  of  ear  protection; 
and  it  is  necessary  to  subtract  the  attenuation  to  be  expected  of  such 
devices  in  determining  the  values  for  S&v  to  be  used  with  these  charts. 

For  convenience.  Table  VI  lists  typical  values  for  such  attenuation  as 
obtained  from  Ref*  2.  It  is  recognized,  of  course,  that  the  actual  reduc¬ 
tion  achieved  for  an  earplug  depends  to  a  great  extent  on  the  fit*  For 
this  reason,  an  ear  muff,  or  combination  of  muff  and  plug,  is  often 
preferred. 

TABLE  VI 

Average  Octave-Band  Noise  Level  Reduction 
in  Bar  Canal  Achieved  by  Use  of  Personal 
Protective  Equipment 


Bands 

(cps) 

Type  of 

600 

1200 

2400 

Equipment 

1200 

2400 

4800 

Headset  Earphone  Covers 

13 

20 

30 

Standard  Ear  Plug  or 

Ear  Muff  alone 

18 

23 

30 

Standard  Ear  Plug  and 

Ear  Muff  together 

28 

37 

40 

2.  The  Noise  Cumulator 

9  14 

Several  years  *-go,  Cox  9  developed  an  instrument  known  as  the 
Noi9e  Cumulator  for  the  specific  purpose  of  studying  noise  exposures 
aboard  aircraft  carriers.  In  the  design  of  this  instrument  che  assumption 
was  made  that  the  total  amount  of  time  a  person  is  exposed  to  noise  of  a 
given  level  and  spectrum  determines  the  seriousness  of  his  exposure 
independent  of  the  time  sequence  of  quiet  and  noisy  periods.  Using  the 
actual  noise  signal  (or  a  tape  recording  of  it)  the  cumulator  was  built 
to  tell  the  total  time,  during  any  period  of  interest,  for  which  the 
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level  in  any  specific  frequency  range  exceeded  each  of  seven  pre- determined 
sound  level  thresholds.  These  thresholds  are  normally  spaced  5  db  apart, 

generally  covering  an  ample  dynamic  range  with  reasonable  definition. 

3  4 

As  previously  discussed,  more  recent  studies  *  show  that  the 
seriousness  of  a  noise  exposure  is  not  completely  independent  of  the  time 
sequence  of  quiet  and  noisy  periods  -  since  TTG  recovery  occurs  during 
quiet  periods.  However,  if  the  changes  in  noise  level  are  fairly  rapid  * 
as  is  characteristic  of  aircraft  take-offs  from  carrier  decks  -  recovery 
during  quiet  periods  becomes  less  important,  and  the  information  supplied 
by  the  cumulator  is  sufficient  for  a  calculation  of  TT^  in  a  manner 
similar  to  that  previously  outlined. 

In  order  to  demonstrate  the  applicability  of  this  instrument  to  pre¬ 
diction  of  TTS^,  we  have  made  use  of  Noise  Cumulator  data  taken  aboard 
U.S.S.  Saratoga  (CVA-60)  at  the  Blast-Shield  Operator's  Station,  Catapult 
No.  1,^  The  authors  of  the  CVA-60  study  identify  this  as  the  noisiest 
position  which  they  monitored. 

Three  frequency  ranges  were  studied:  full  band-width,  300-1200  cps, 
and  1200-20,000  cps.  Data  was  recorded  on  eight  magnetic  tapes  ("sub¬ 
samples"),  each  approximately  50  minutes  long.  Five  of  these  were  defined 
as  "loud",  i.e.,  the  overall  levels  exceeded  130  db  more  than  2%  of  the 
time.  Our  analysis  is  based  on  these  sub-samples,  under  the  assumption 
that  these  samples  represent  continuously  noisy  conditions  without  appreci¬ 
able  quiet  periods  for  recovery. 

Figure  10  shows  the  percentage  of  time  that  the  sound  level  spent 
in  each  5  db  increment  for  each  of  the  frequency  ranges  monitored.  It  is 
immediately  clear  that  most  of  the  high-level  sound  energy  was  in  the  1.2- 
20  KC  range;  this  is  shown  by  the  similarity  of  the  ful 1- bandwidth  and  the 
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1.2-20  KC  plots.  The  distribution  curve  for  300-1200  cps  peaked  some 
15  to  20  db  lower. 

This  would  not  be  expected  if  the  noise  measured  was  due  chiefly  to 
Jet  engine  exhausts.  The  predominance  of  high  frequency  energy  may  be  due 
to  compressor  noise  from  aircraft  awaiting  takeoff,  and  to  noise  from  the 
catapult  itself. 

As  discussed  previously,  TTS^  at  4000  cps  is  chiefly  influenced  by 
noise  in  the  1200-2400  cps  and  the  2400-4800  cps  octave  bands.  TTS^  at 
2000  cps  is  governed  by  the  600-1200  and  1200-2400  cps  bands.  The  CVA-60 
data  does  not  give  this  information,  but  we  can  make  some  assumptions  which 
will  allow  us  to  estimate  how  the  measured  energy  should  be  divided  up. 

The  fact  that  this  involves  some  educated  guessing  points  out  the  desira¬ 
bility  of  using  filters  tailored  to  the  three  octaves  influencing  TTS  at 
the  two  frequencies  of  most  interest. 

Let  us  assume  that  the  1.2-20  KC  band  noise  is  mainly  due  to  engine 
compressor  noise.  As  such,  it  may  be  expected  to  have  approximately  equal 
energy  in  the  1200-2400  and  2400-4800  cps  bands,  the  energy  content  being 
a  function  of  engine  type  and  speed.  Energy  at  higher  frequencies  will  be 
negligible  compared  to  these  bands. 

We  further  assume  that  the  noise  in  the  300-1200  cps  band  is  due  to 
turbojet  engine  exhausts  (at  different  power  settings),  plus  miscellaneous 
sources  such  as  tractors  and  the  public  address  system.  We  can  assume  that 
the  levels  in  the  300-600  and  600-1200  cps  octaves  which  make  up  this  band 
are  about  the  same. 
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With  these  assumptions  we  can  say  that  the  600-1200  cps  octave  level 
is  3  db  below  the  300-1200  cps  band  level,  and  that  the  1200-2400  and 
2400-4800  cps  levels  are  each  3  db  below  the  1*2-20  KC  band  level.  Taking 
the  average  level  in  each  5  db  increment  as  being  at  about  the  middle  of 
that  range,  we  can  now  calculate  TTS  for  a  subject  wearing  V-51R  ear  plugs, 
having  attenuations  of  18  db  in  the  600-1200  cps  band,  23  db  at  1200-2400 
cps  and  30  db  at  2400-4800  cps.  We  make  use  of  the  growth  equations  in 
Table  II  with  5  db  corrections  for  pure  tones  in  the  1200-2400  and  2400- 
4800  cps  bands. 

Figure  11  details  the  procedure  in  calculating  TTS^  at  4  KC  due  to 
the  1200-2400  cps  band,  namely  22.6  db  for  50  minutes  of  exposure  during 
noisy  periods.  A  similar  calculation  employing  the  2400-4800  cps  octave 
gives  a  TTS^  of  20,1  db,  the  total  TTS^  at  4  KC  being  the  sum  of  these, 
or  24.6  db.  The  TTS^  at  2000  Cj.  calculated  in  a  similar  manner,  totals 
9.3  db.  These  values  of  TTS^  corroborate  the  conclusion  of  reference  15 
as  to  the  desirability  of  using  both  insert  and  muff-type  hearing  protec¬ 
tion  simultaneously. 
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(Add  5  db  for  pure  tones) 
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(Add  5  db  for  pure  tones) 
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OCTAVE-BAND  SOUND  LEVEL  RANGE  IN  DB 
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3.  Speech  Communication  in  the  Presence  of  Noise 

In  order  to  estimate  the  effects  of  noise  on  speech  communication, 
the  level  and  spectrum  of  the  noise  are  often  described  in  terms  of  a 
speech- inter ference-level^*^  By  definition,  the  speech- inter ference- 
level  is  the  arithmetic  average  sound  level  (at  the  point  of  interest) 
of  several  specified  octave-bands.  Tabulations  are  then  given  of  the 
maximum  distance  for  reliable  conversation  as  a  function  of  speech-in¬ 
ter  ference- level  (SIL)  and  voice  level-the  latter  being  specified  simply 
as  normal,  raised,  etc.  There  is  some  disagreement,  however,  as  to 
which  octave  bands  should  be  employed  in  calculating  SIL.  For  example, 
Beranek16,17  uses  the  3  bands  600-1200,  1200-2400,  and  2400-4800  cps ; 

Scrasberg^®  suggests  the  4  bands  300-600,  600-1200,  1200-2400,  and 
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2400-480C  ops;  while  recent  studies  of  Klumpp  and  Webster  seem  to 
correlate  better  with  the  3  bands  300-600,  600-1200,  1200-2400  cps. 

Also,  slight  modifications  of  the  SIL  criteria  should  be  made  for  the 
several  calculation  methods. 

Actually,  the  concept  of  speech- interference-level  is  based  upon 

a  simplification  of  a  method  originally  advanced  by  French  and  Stein- 
20-22 

berg  for  the  calculation  of  Articulation  Index  (AI>  from  the  speech- 

to-noise  ratio  in  20  frequency  bands.  These  bands  are  intended  to  be 

of  equal  importance  from  the  standpoint  of  their  contribution  to  speech 
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intelligibility.  Recent  tests  confirm  that  the  AI  calculation  is 

indeed  a  good  method  of  predicting  speech  intelligibility. 

Since  AI  calculation  is  adequate,  the  very  fact  that  disagreement 
has  arisen  over  which  octaves  should  be  employed  in  the  simplified  SIL 
approach  indicates  that  we  are  attempting  to  oversimplify  the  problem. 
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This  Is  undoubtedly  due  partially  to  the  fact  that  the  octaves  employed 
do  not  contain  all  of  the  original  20  bands;  nor  do  the  various  octave- 
bands  contain  equal  numbers  of  the  20  bands  of  equal  importance.  Another 
factor  which  is  usually  not  considered  at  all  is  the  fact  that  SIL 
criteria  as  normally  used  correspond  to  an  AI  of  approximately  0.40-which 
although  sufficient  for  sentence  intelligibility  is  not  always  enough 
for  high  intelligibility  for  isolated  words.  (There  is  one  paper  in  the 
literature  which  does  attempt  to  extend  an  SIL  calculation  to  other 
values  of  AI,  )  In  addition,  no  SIL  approach  to  speech  intelligibility 
has  been  proposed  chat  takes  into  consideration  many  other  factors  such 
as  reverberation,  frequency  and  amplitude  distortion  of  speech  by  a  PA 
system,  intermittent  noise,  spread  of  masking,  vocal  effort,  ear  over¬ 
load.  j,  etc, -all  of  which  may  be  brought  into  the  picture  when  AI  cal¬ 
culation  is  employed.  Hence,  from  a  technical  standpoint,  it  seems  that 
the  uae  of  AI  should  be  favored  over  SIL.  Of  course,  the  primary  reason 
why  SIL  has  gained  favor  is  its  basic  simplicity, 

A  recent  paper  by  Kryter^/,iCH  seems  to  provide  an  answer  to  this 
dilemma.  He  proposes  the  use  of  the  six  octave  bands  150-300  cps  thru 
4800-9600  cps  with  a  weighting  factor  for  each  band  that  depends  upon 
how  many  of  the  original  20  bands  fall  within  the  octave.  Taking  these 
factors  as  n/600  (when  n  is  the  number  of  speech  bands  vhereln  the  octave) 
we  may  quite  readily  compute  AI  for  a  given  speech- to-noise  ratio  by 
adding  the  contributions  of  the  6  octaves.  Table  VII  lists  these 
weighting  factors  and  shows  the  basic  method  of  AI  calculation.  In 
such  calculations,  the  maximum  effective  speech- to-noise  ratio  is  taken 
as  30  db,  while  ratios  of  0  db  or  less  contribute  nothing. 
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Table  VII 


Octave  band 

150 

300 

600 

1200 

2400 

4800 

Cps 

300 

600 

1200 

2400 

4800 

9600 

1*  Weighting  Factor 

.0013 

.0042 

.0067 

.0105 

,0089 

0017 

2.  Speech  Peak  Level 

77 

80 

75 

70 

64 

58 

3.  Noise  Level 

69 

67 

62 

53 

51 

48 

4*  S/N  ratio,  db  (#2-#3) 

8 

13 

13 

17 

13 

20 

5.  AI  contribution  (#1  x  #4) 

.0104 

.0546 

.0871 

.1786 

.1156 

.0340 

6.  AI  -  sum  of  line  #5  z  0.48 

From  Figure  12,  this  corresponds  to  about  97%  sentence  intelligibility 
28  20 

,  Figure  12  also  shows  average  intelligibility  scores  for  various 
lists  of  words  and  syllables*  Such  calculations  agree  well  with  the  more 
involved  20  band  procedures . ^ 

Thus,  a  knowledge  of  the  speech- to-noise  ratio  in  6  octaves  provides 

a  better  picture  of  speech  intelligibility  than  can  be  obtained  from  SIL 

alone*  Actually,  the  calculation  is  not  much  more  difficult;  noise  levels 

are  needed  in  3  or  4  of  these  6  octaves  in  order  to  compute  SIL,  and 

tests  have  shown  that  a  generalized  speech  spectrum  is  sufficient.  It 

Is  only  necessary  to  determine  the  overall  speech  level-which,  of  course, 

varies  with  voice  level,  distance,  and  use  of  PA  system,  if  any,  Table 

VIII  gives  typical  levels  for  speech  pea ks  (male  voices)  when  the  long 
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term  rms  overall  speech  level  is  65  db.  9  9  This  is  approximately 

a  normal  speech  level  at  a  distance  of  1  meter.  In  practice,  the  long 
term  rme  level  can  be  estimated  by  subtracting  3  db  from  the  arithmetic 
average  of  the  peak  readings  of  a  sound  level  meter  set  for  slew  response 
on  the  "C"  scale.  Actual  speech  peaks  are  normally  about  12  db  above 
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the  long  term  rms  level-or  9  db  above  the  readings  of  such  a  meter.  If 
a  sound  level  meter  la  not  available,  one  may  make  a  rough  estimate  of 
speech  level  from  the  fact  that  a  "raised  voice"  is  about  6  db  louder 
than  a  "normal"  voice,  and  a  "very  loud"  voice  about  6  db  louder  than 
a  "raised"  voice.  Such  assumptions  are  exactly  those  used  in  the  sim¬ 
plified  SIL  approach.^  Obviously,  of  course,  speech  level  will  decrease 
with  increased  distance-6  db  per  doubling  of  distance  in  a  free  field 
and  somewhat  less  than  that  in  a  normal  room.  (The  SIL  approach  simply 
uses  the  6  db  per  doubling  figure.) 

Table  VIII 


Octave  Band 

150 

300 

600 

1200 

2400 

4800 

Cps 

300 

600 

1200 

2400 

4800 

9600 

Typical  Speech  Peak  Levels 
(65  db  rms  overall) 

71 

74 

69 

64 

58 

52 

When  PA  systems  are  used,  the  long  term  average  rma  speech  level 
must  be  measured  at  the  location  of  interest.  Table  VIII  may  then  be 
used  to  determine  the  speech  spectrum.  For  example,  if  the  average  rms 
level  were  found  to  be  87  db,  the  speech  spectrum  would  be  estimated  by 
adding  87-65*  22  db  to  the  levels  of  Table  VIII. 

The  previous  discussion  is  based  on  the  assumption  that  the  calker 
and  listener  cannot  see  each  other,  that  reverberation  time  is  zero, 
that  the  noise  is  steady-state,  that  the  communications  link  is  distor¬ 
tionless,  etc.  In  practice,  all  such  requirements  are  not  likely  to  be 
fulfilled  in  a  given  situation.  Kryter  has  discussed  proposed  modifiers 
to  the  basic  technique  to  permit  extension  to  these  situations.  These 
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will  be  briefly  discussed.  All  of  these  factors  are  normally  neglected 
in  the  SIL  approach. 

Figure  13  shows  the  basic  work-sheet  for  calculation  of  Al  by  octave 
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bands,  including  the  spectrum  for  normal  male  speech.  if  the  speech 
is  electronically  amplified,  there  is  a  maximum  tolerable  sound  level 
which  should  not  be  exceeded  since  it  overloads  the  hearing  mechanism 
and  does  not  contribute  to  intelligibility.  If  the  speech  link  does 
not  have  flat  frequency  response,  the  effective  speeoh  spectrum  should 
be  adjusted  accordingly.  Another  factor  which  should  be  considered  for 
loudspeaker  presentation  is  indicated  in  Table  IX.  This  is  a  correction 
for  degradation  in  speech  intelligibility  at  high  levels  in  a  partially 
reverberent  room.27,30  ic  should  not  be  applied  outdoors. 

Table  IX 


Overall  rms  speech 
level 

Amount  to  be  subtracted 
from  speech  level 

Resulting  Effec¬ 
tive  speech  level 

85  db 

0  db 

85  db 

90 

2 

88 

95 

4 

91 

100 

7 

93 

105 

11 

94 

no 

15 

95 

115 

19 

96 

120 

23 

97 

125 

27 

98 

130 

30 

100 

Thus,  little  is  to  be  gained  by  increasing  indoor  speech  levels  above 
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about  100  dbf 


Speech  communications  systems  sometimes  employ  peak  clipping  to 
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decrease  the  peak-to-rms  ratio,  thereby  improving  AI.  Figure  14 

shows  the  effective  increase  in  speech  level  to  be  added  to  the  speech 

spectrum  normally  used  with  Fig.  13  -  as  a  function  of  the  amount  of 
31  27 

clipping.  *  Twenty- four  db  of  clipping  is  equivalent  to  12  db 
greater  speech  level. 

If  the  talker  is  not  speaking  in  a  reasonably  normal  voice,  a 
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further  adjustment  of  the  AI  Is  necessary.  If  a  very  low  or  very 
loud  tone  is  used,  the  effective  speech  level  -  to  be  used  in  conjunc¬ 
tion  with  Fig.  13  -  must  be  lower  than  the  true  speech  level  since 
vocal  quality  will  differ  from  normal  speech.  Figure  15  shows  this 
relationship  *  ^ *  27 

Studies  also  show  that  very  high  background  noise  levels  provide 

more  speech  masking  for  a  given  speech- to-noise  ratio  than  low  back- 

ground  noises  with  the  same  spectrum.  ^  This  introduces  a  correction 

to  the  background  noise  for  each  frequency  band  which  is  determined  by 

the  excess  of  the  background  noise  above  the  threshold  of  audibility 

In  that  band.  Table  X  and  Figure  13  both  show  the  average  threshold 

of  audibility,  while  Table  XI  lists  the  correction  to  be  applied  to 
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the  noise  spectrum.  * 


Table  X 


Octave -band 

150 

300 

600 

1200 

2400 

4800 

Ops 

300 

600 

1200 

2400 

4800 

9600 

Threshold  of  audibility 

26 

17 

12 

9 

-S 

13 
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Table  XI 


Band  Level  minus  Correction  to  be  added  to 

Threshold  Level  Noise  Level  for  this  Band 


80  db 
85 
90 
95 
100 
105 

no 

115 

120 

125 

130 

135 

140 

145 

150 


0  db 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 


If  the  noise  is  predominately  of  a  narrow  bandwidth  character 

still  another  correction  should  be  applied  for  upward  and  downward 

spread  of  masking  beyond  the  actual  bandwidth  of  the  noise  per  se. 

35,27  However,  for  broadband  noises  -  as  are  characteristic  of  jet 
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aircraft,  etc.  -  such  corrections  are  much  less  important  and  so 
will  not  be  considered  in  detail  in  this  report. 

If  the  background  noise  is  interrupted,  the  calculated  AI  is 
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As  would  be  expected, 


adjusted  by  means  of  Figs.  16  and  17. 
interruption  of  the  noise  increases  AI;  the  amount  of  increase  is 
strongly  influenced  by  interruption  rate.  Aircraft  take-off  from  a 
carrier  deck  is  normally  at  too  low  a  frequency  to  Increase  AI  appreci¬ 
ably.  Apparently  a  10  cps  interruption  rate  is  optimum  to  permit  the 
listener  to  hear  enough  of  each  word  unmasked  to  piece  together  the 
whole.  At  very  high  interruption  rates  the  noise  acts  like  a  contin¬ 
uous  noise  of  somewhat  lower  level. 

Reverberation  acts  to  decrease  AI  as  shown  by  Fig.  18. ^7>27  ^ 

a  carrier  flight-deck  the  reverberation  time  would  be  aero,  but  would 
not  be  negligible  in  an  enclosed  area. 

In  the  same  way  that  a  deaf  person  can  lip-read,  visual  cues 

give  an  improvement  in  AI  when  the  listener  can  see  the  talker.  As 

would  be  expected,  the  poorer  the  AI  without  visual  cues,  the  more 
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they  help,  as  in  Fig,  19*  * 

The  four  examples  given  below  illustrate  the  calculation  and  use 
of  AI  for  hypothetical  situations.  They  should  also  serve  to  show  the 
superiority  of  an  AI  calculation  over  the  normal  SIL  procedures. 

Example  Is  The  noise  levels  at  a  particular  point  of  interest  due  to 
a  relatively  distant  group  of  aircraft  (perhaps  at  idle  power)  is  as 
tabulated  in  line  1  of  Table  XII.  A  nearby  loudspeaker,  operating  from 
a  PA  amplifier  etc.,  delivers  an  average  speech  level  -  as  read  by  the 
peaks  of  the  slow  response  of  a  sound  level  meter  on  the  C  weighting 
network  *  of  123  db.  The  PA  system  has  a  flat  frequency  response  (in¬ 
cluding  microphone)  except  at  high  frequencies  -  it  being  2  db  down  in 
the  2400-4800  cps  range  and  5  db  down  for  4800-9600  cps. 
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As  previously  discussed,  the  long  term  average  rms  speech  level 
may  be  estimated  by  subtracting  3  db  from  the  sound  level  meter  read* 
ings  obtained  as  described  above.  The  result  is  120  db.  Table  VIII 
gives  a  typical  peak  speech  spectrum  for  a  65  db  rms  level.  Hence, 
we  need  only  add  120-65  3  55  db  to  these  figures  to  obtain  the  peak 
speech  spectrum  for  a  flat  PA  system.  The  result  is  given  in  line  2 
of  Table  XII.  Applying  the  minor  high  frequency  less  corrections  for 
the  PA  system,  we  obtain  line  3  of  this  tabulation.  (If  such  a  pro¬ 
cedure  is  sufficient  to  change  the  overall  level  appreciably,  all 


numbers  should  be  increased 

quired  to  maintain  overall 

.  or  decreased  by  tbe 

level  ,) 

Table  XII 

constant 

amount 

re- 

Octave  Bands 

150 

300 

600 

1200 

2400 

4800 

Cps 

300 

600 

1200 

2400 

4800 

9600 

1. 

Meas.  Noise  Level 

103 

109 

111 

113 

106 

100 

2. 

Speech  Peaks  (uncorrected 
for  PA  system) 

126 

129 

124 

119 

113 

107 

3. 

Speech  Peaks  (corrected 
for  PA  system) 

126 

129 

124 

119 

111 

102 

4. 

Audibility  Threshold 

26 

17 

12 

9 

3 

13 

5. 

Line  1  -  Line  4 

77 

92 

99 

104 

103 

87 

6. 

Correction  to  Noise 

- 

2 

4 

5 

5 

- 

7. 

Corr.  Noise  Spectrum 

103 

111 

115 

118 

111 

100 

8. 

Line  3  -  Line  7 

23 

18 

9 

1 

0 

2 

9. 

Weighting  Factor 

.0013 

0042 

.0067 

.0105 

.0089 

.0017 

10. 

11. 

Contributions  to  AI 

A I  3  sum  of  Line  10  o  0.18 

.030 

-37 

.076 

.060 

.011 

- 

.003 

-37- 


Since  Line  3  does  not  exceed  the  maximum  tolerable  levels  for  unclipped 
speech  -  as  shown  on  Figure  13  -  no  further  speech  corrections  are  re¬ 
quired. 

Line  4  gives  the  audibility  threshold  as  in  Table  X,  and  the  in¬ 
dicated  corrections  for  high  level  masking  (line  6)  are  obtained  from 
Table  XI. 

Subtracting  line  7  from  line  3,  we  obtain  the  effective  speed- 
to-noise  ratio  for  each  band.  These  must  be  multiplied  by  the  weighting 
factors  (line  9)  from  Table  VII  to  give  the  contributions  of  the  in¬ 
dividual  octave  bands  to  the  AI,  The  total  AI  is  the  sum  of  these  con¬ 
tributions,  or  0,18. 

From  Fig.  12,  this  would  be  expected  to  provide  only  about  757* 
sentence  intelligibility,  and  very  poor  intelligibility  for  individual 
words . 

Example  2:  The  noise  Level  is  the  same  as  for  example  1.  However, 

12  db  speech  clipping  is  employed  in  the  PA  system  with  12  db  post¬ 
clipping  gain  to  achieve  the  previous  peak  amplitudes.  . 

Referring  to  Fig.  14,  the  effective  speech  peaks  have  been  in¬ 
creased  9  db  due  to  clipping  and  post-clipping  gain.  Hence,  the 

i 

speech  levels  are  as  indicated  in  the  first  line  of  Table  XIII  -  being 
obtained  by  adding  9  db  to  the  corresponding  figures  of  line  3,  Table 
XII.  In  the  600-1200  band,  these  effective  speech  peak  levels  exceed 
the  maximum  tolerable  levels  of  Fig.  13  for  12  db  clipping  by  a  small 
amount.  Hence,  the  speech  peaks  are  corrected  as  indicated  in  line  2. 
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Table  XIII 


Octave  Bands 

130 

300 

600 

1200 

2400 

4800 

Cp8 

300 

600 

1200 

2400 

4800 

9600 

1. 

Effective  Peaks 

135 

138 

133 

128 

120 

111 

2. 

Eff.  Speech  Peaks  corrected 
for  maximum  level 

135 

138 

132 

128 

120 

111 

3. 

Corr.  Noise  Spectrum 

103 

111 

115 

118 

111 

100 

4, 

Line  2  -  Line  3 

32 

27 

17 

10 

9 

11 

5. 

Eff.  Speech- to-noise 

30 

27 

17 

10 

9 

11 

6. 

Weighting  Factor 

0013 

0042 

.0067 

.0105 

.0089 

.0017 

7. 

Contribution  to  AI 

039 

113 

.114 

.105 

.080 

.019 

8. 

AI  s  sum  of  Line  7  *  0.47 

Line  3  shows  the  corrected  noise 

spec 

crum,  as 

taken 

from  line 

7 

of  Table  XII.  The  calculation  than  proceeds  as  before,  except  chat  in 
the  150-300  cps  band  it  ts  necessary  to  make  use  of  the  fact  that  the 
effective  speech- to-noise  ratio  should  be  limited  to  30  db. 

The  AI  is  now  computed  as  0*47,  which  from  Table  12  corresponds 
to  977o  sentence  intelligibility  and  947«  word  intelligibility  for  the 
250  word  list.  Thus,  a  worthwhile  improvement  has  been  obtained  by 
speech  clipping. 

Example  3:  Assume  the  locale  to  be  a  conference  room  at  some  distance 
from  flight  deck  operations.  Assume  the  noise  levels  are  as  given  in 
line  1  of  Table  XIV,  and  that  the  room  has  a  reverberation  time  of 
about  1  second  (at  500  cpa)  Further  assume  that  no  PA  system  Is  being 
used,  but  that  the  speaker  is  talking  in  a  rather  loud  voice.  For 
example,  take  the  long  term  rms  speech  level  as  80  db  at  1  meter- though 
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only  70  db  at  a  typical  location  of  interest.  (The  sound  level  meter 
readings  on  C  -  slow  would  be  about  3  db  higher  for  peaks)  It  may  also 
be  assumed  that  the  audience  is  watching  the  speaker,  and  hence  is  aided 
somewhat  by  visual  cues. 

From  Fig.  15  it  is  seen  that  high  vocal  effort  is  limiting  the 
speech  intelligibility  slightly,  and  that  to  correct  for  this  a  2  db 
subtraction  should  be  made  from  the  actual  speech  levels.  Hence,  the 
effective  speech  level  at  a  typical  position  of  interest  is  70  -  2  ®  68 
db.  Then  from  VIII,  the  speech  peak  spectrum  is  as  Indicated  In  line  2 
of  Table  XIV. 


Table  XIV 


Octave  Bands 

150 

300 

600 

1200 

2400 

4800 

Cpa 

300 

600 

1200 

2400 

4800 

9600 

1. 

Noise  Spectrum 

67 

57 

52 

49 

43 

40 

2. 

Corrected  Speech  Peaks 

74 

77 

72 

67 

61 

55 

3. 

Line  2  -  Line  1 

7 

20 

20 

18 

18 

15 

4. 

Weighting  Factor 

.0013 

.0042 

.0067 

.0105 

.0089 

.0017 

5. 

Contribution  to  AI 

.009 

.084 

.134 

.189 

.160 

.026 

6.  AI  =  sum  of  Line  5  ■  0.60 

7.  AI  corrected  for  1  sec  reverberation  time  -  0.50 

8.  Effective  AI  with  visual  ever*  -  0.63 

Since  neither  noise  nor  speech  levels  are  excessive,  no  further 
corrections  are  necessary,  and  the  calculation  proceeds  in  a  normal 
fashion  as  outlined  in  Table  XIV,  thru  line  6. 

Now,  however,  a  correction  is  applied  for  the  1  second  assumed 
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reverberation  time.  From  Fig.  18,  this  is  expected  to  reduce  AI  by  0.1. 
Fig.  19  indicates  that  this  reduction  is  more  than  compensated  for  by 
the  effect  of  visual  cues  -  the  final  result  being  AI  =  0.63. 

Example  4:  Assume  the  same  room  as  in  example  3,  except  that  the  noise 
spectrum  has  been  raised  25  db  and  a  flat  response  PA  system  used  to 
increase  the  speech  level  25  db  without  clipping  -  thus  maintaining  the 
same  speech-to-noise  ratio. 

The  speech  level  is  now  70  /  25  *  95  db  rms.  Hence,  from  Table  IX 

the  effective  level  (in  a  room)  is  but  91  db.  The  speech  peak  spectrum 

is,  therefore,  91  -  65  =  26  db  above  that  of  Table  VIII.  This  is  shown 

in  line  1  of  Table  XV.  Line  2  shows  the  noise  spectrum  as  described 
above.  Since  neither  speech  nor  noise  is  excessive  (See  Fig.  13  and 
Tables  X  and  XI)  the  calculation  now  proceeds  in  a  normal  manner  as  in¬ 
dicated  In  Table  XV.  The  results  indicate  only  a  slight  drop  in  in¬ 
telligibility  from  example  3  to  example  4.  However,  if  in  example  3 
the  speaker  had  been  talking  with  a  normal  vocal  effort,  the  decrease 
In  AI  due  to  raising  both  speech  and  noise  by  equal  amounts  would  have 
been  somewhat  greater. 


Octave  band 
Cps 

1.  Effective  Speech  Peaks 

2.  Noise  Spectrum 

3.  Line  1  *  Line  2 

4.  Weighting  Factor 


Table  XV 


150 

300 

300 

600 

97 

100 

92 

82 

5 

18 

0013 

.  0042 
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600 

1200 

1200 

2400 

95 

90 

77 

74 

18 

16 

0067 

.0105 

2400 

4800 

4800 

9600 

84 

78 

68 

65 

16 

13 

0089 

.0017 

Table  XV  con1 t 


5.  Contributions  to  AI  .007  ,076  .121  .168  .142 

6.  AI  =  sum  of  line  5  c  0.54 

7.  AI  corrected  for  1  sec.  reverberation  time  :  0.44 

8.  Effective  AI  with  visual  cues  •  0.60 


.022 
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SYLLABLES,  WORDS  OR  SENTENCES  UNDERSTOOD  -  PERCENT 
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DB  INCREASE  IN  LONG  TERM  RMS  PRESSURE  OF  PEAK-CLIPPED 
SPEECH  RELATIVE  TO  LONG  TERM  RMS  PRESSURE  OF  UNCLIPPED 
S PERCH  OF  EQUAL  PEAK  AMPLITUDE 


AMOUNT  OF  PEAK  CUPPING  and  post-clipping  in  db  (re  0.001 

PROBABILITY) 


FIGURE  14 
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ACTUAL  SPEECH  LEVEL  1  METER  FROM  TALKER 
LONG  TERM  RMS  IN  DB  RE  0.0002  MtCROBAR 


ORRECTION  TO  BE  ADOEO  TO  AI 
OMPUTED  ON  BASIS  OF  STEADY 
TATE  NOISE 


oow 


FIGURE 


THE  ORDINATE  SHOWS  A  CORRECTION  TO 
BE  APPLIED  TO  THE  ARTICULATION  INDEX 
COMPUTED  ON  THE  ASSUMPTION  THAT  A 
MASKING  NOISE  IS  STEADY  STATE  FOR 
VARIOUS  NOISE -TIME  FRACTIONS.  (4) 

THE  CORRECTED  AI  CAN  NOT  EXCEED 
ONE  (I.O) 
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EFFECTIVE  AI 


FIGURE  17 
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AMOUNT  TO  BE  SUBTRACTED  FROM  AI 


0  Ig.  1  -  1  -  1  ■ .  1 - 1  -  1  -  -  ■  1 

01  23456789 

REVERBERATION  TIME  IN  SECONDS 


FIGURE  18 
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FIGURE  19 
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IV,  EFFECTS  OF  NOISE  ON  STRUCTURES 


1.  Introduction 

The  high  noise  levels  produced  by  jet  and  rocket  engines  has  led  to 
fatigue  failure  of  structural  parts  and  to  excessive  vibration  resulting 
in  malfunction  of  equipment.  Great  effort  has  been  made  by  many  people 
to  analyze  and  find  solutions  in  various  phases  of  the  problem.  In  spite 
of  the  vast  literature  on  the  subject  there  is  still  no  reliable,  widely 
used  method  for  specifying  the  allowable  noise  levels  for  a  given  aircraft 
or  conversely  to  define  a  structure  to  withstand  a  given  noise  level. 

Very  considerable  progress  has  been  made  in  various  areas  but  the 
problem  is  so  complex  that  much  remains  to  be  done.  There  is  a  need  for 
a  more  accurate  method  for  calculating  the  noise  pressure  levels  in  the 
near  field  of  a  jet.  The  space  correlations  of  the  pressures  which  are 
important  in  determining  the  coupling  between  the  pressure  and  the 
structure  require  more  study.  Although  methods  for  calculating  the 
resonant  frequencies  of  many  types  of  structures  have  been  used  with  some 
success,  the  computation  of  the  response  and  the  stress  or  strain  in 
practical  structures  is  quite  difficult  and  unreliable.  The  damping  is 
of  paramount  Importance  in  determining  the  response  at  resonance  and  is 
usually  measured  or  estimated  rather  than  calculated.  The  mechanism  of 
fatigue  of  structures  is  not  well  understood  and  there  is  considerable 
scatter  in  fatigue  data  and  even  some  contradictory  results. 

Fatigue  tests  on  complete  structures  using  jet  noise  sources  are 
quite  expensive  and  time  consuming.  Therefore,  no  large  amounts  of  data 
have  been  made  available  by  such  means.  Instead,  most  fatigue  tests  are 
conducted  on  small  structures,  test  panels  or  even  small  beams  with  loud 
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speakers  and  sirens  or  mechanical  load  cycling  machines  are  used  to  apply 
the  stress  loads.  The  data  must  be  correlated  with  prototype  data  under 
real  service  conditions.  Fitch  and  his  associates  have  reported  (Ref.  39) 
what  appears  to  be  a  very  promising  approach  to  the  problem.  They  combine 
engineering  analysis  with  accelerated  life  testing  of  a  few  key  structures 
to  determine  whether  the  structures  have  adequate  design  margins.  This 
method  outlined  below  will  be  discussed  more  fully  in  later  paragraphs. 

The  analysis  of  the  problem  includes  the  following  steps: 

(a)  Determine  the  acoustic  environment  under  various  operating 
conditions*  Measured  data  of  course  are  preferred  but,  if 
necessary,  predictions  can  be  made  from  engine  parameters. 

(b)  Analyze  the  service  usage  and  missions  of  the  aircraft  to 
determine  the  proportion  of  the  design  life  th zz  the  aircraft 
is  exposed  to  each  of  the  important  acoustic  loads. 

(c)  Build  test  specimen  to  simulate  each  critical  structure  and 
test  to  failure  using  sinusoidal  excitation. 

(d)  Determine  the  relation  between  statistical  response  to  random 
excitation  and  sinusoidal  excitation  for  simple  structures. 

(e)  Using  the  data  from  the  test  specimen  and  a  random  S-N  curve, 
find  the  allowable  equivalent  random  sound  pressure  level  which 
would  give  the  required  design  life  for  the  critical  operations 
determined  in  the  service  usage  analysis  of  step  (b) . 

(f)  Compare  the  allowable  sound  pressure  level  with  the  actual  sound 
pressure  to  see  if  adequate  margins  exist. 

Belcher  (Ref.  40)  and  McGowan  (Ref.  41)  have  reported  the  development 
of  design  charts  suitable  for  use  by  the  average  structural  design  engineer. 

The  charts  in  effect  are  nomograms  relating  sound  pressure,  fatigue  life, 

-52- 


resonant  frequency  and  the  dimensions  of  several  common  panel  constructions. 


2.  Determination  of  Sound  Levels 

The  most  severe  acoustic  environment  for  the  aircraft  structures 
probably  occurs  at  maximum  power  during  static  ground  operation.  As  the 
vehicle  gains  speed  and  altitude  the  engine  noise  decreases.  Although 
boundary  layer  noise  Increases  with  speed  it  is  seldom  as  important  as  the 
maximum  engine  noise. 

Measurements  of  the  overall  sound  pressure  and  octave  band  frequency 
analyses  near  the  aircraft  are  the  preferred  method  of  describing  the 
acoustic  environment  but  it  can  be  estimated  using  the  methods  given  by 
Fitch  (Ref.  39). 

The  procedure  is  a  straight  forward  step  by  step  calculation  using 
standard  reference  contours  and  correction  factors  based  on  empirical  data. 
The  corrections  are  functions  of  the  exit  diameter  and  expanded  exhaust 
velocity. 

The  steps  in  the  calculation  are  outlined  below  and  illustrated  at 
Military  Rated  Thrust  (MRT)  and  Maximum  Afterburner  (Max  A/B)  with  data 
from  Table  XVI  (Ref.  13,  Appendix  D) . 

1.  Calculate  the  effective  exhaust  velocity  of  the  engine. 

V  ■  tg  t  *  Thrust/engine*  pounds 

6  g  *  32.2  ft. /sec. 

w  =  weight  flow  rate,  pounds/sec. 

MRT:  V  =  8310  x  32.2  ■*  1871  ft. /sec. 

6  143 

Max  A/B:  V  -  12630  x  32.2  =  2745  ft. /sec. 
e  143 
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TABLE  XVI 


J79-2  Engine  Operating  Data 
Installed  in  A3J-1  Aircraft,  No  Compressor  Bleed 
(From  Appendix  D,  Reference  13) 


Idle 

M.R.T. 

Max  A/B 

Thrust  per  Engine,  lb. 

250 

8310 

12630 

Weight  Flow,  lb. /sec. 

47.5 

143 

143 

2 

Nozzle  Area,  in. 

585 

350 

604 

Nozzle  Diameter,  ft. 

2.275 

1.760 

2.310 

Nozzle  Pressure  Ratio 

1.043 

2.170 

1.935 

Nozzle  Total  Temp.,  °R 

959 

1580 

3500 

o 

Nozzle  Static  Temp.,  R 

950 

1360 

3040 

Nozzle  Gas  Velocity,  ft. /sec. 

170 

1871 

2745 

Ratio  of  Specific  Heats 

1.40 

1.33 

1.25 

2.  Calculate  and  add  the  change  in  sound  pressure  level  to  each 

reference  contour,  utilizing  the  effective  velocity  and  the 

velocity  exponent,  n,  for  each  contour,  Fig.  20. 

ASPL  «  10  n  log  \ 

1850 

MRT:  ASPL  =  10  n  log  1871  =  0.02  n  (negligible) 

1850 

Max  A/B:  ASPL  =  10  n  log  2745  =  1.718  n 

1850 

n  ASPL 

4  7 

5  9 

6  10 

7  12  These  corrections  are  added  to  Fig.  20. 
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3.  Multiply  the  dimensionless  parameters  x/d  and  y/d  in  Fig.  20 
by  the  exhaust  exit  diameter  to  adjust  sound  pressure  level 
contours  to  the  air  vehicled  dimensions. 

For  both  MRT  and  Max  A/B  the  exit  diameter  is  1.76  ft.  hence 
x  *  1.76  |xj  and  y  *  1.76  |  . 

4.  Shift  the  contours  in  Item  3  downstream  a  distance  when 

extrapolating  to  supersonic  exhaust  velocities. 

Ax  =  6.5  D  (M  -  l)^ 
e  '  e 

M  is  very  nearly  1.0  for  both  MRT  and  Max  A/B  and  the  correction 
e 

is  negligible. 

5.  Rotate  the  contours  in  Item  4  through  the  angle  A 0  which  is 
determined  from  Fig.  21  about  the  point  on  the  jet  axis  Ax 
downstream. 

MRT:  At  V  «  1871,  A0  «  0 
e 

Max  A/B:  At  Vg  =  2745,  A0  =  14° 

6.  Calculate  the  frequency  spectra  for  any  position  in  the  near 
sound  field,  from  Fig.  22,  from  knowledge  of  the  over-all  free- 
field  sound  presrure  level  in  Item  5  and  the  jet  velocity  and  exit 
diameter. 

The  first  five  steps  were  carried  out  as  outlined  above  for  the  A3J 

at  MRT  and  Max  A/B.  The  sound  pressure  levels  are  free- field  values  for 

one  engine  operating  above  aground  plane.  Since  the  A3J  has  two  engines 

3db  was  added  to  all  the  contours.  Another  3db  was  added  where  the  contours 

were  superposed  on  the  aircraft  outline  to  account  for  reflections  from 

the  aircraft.  The  resulting  contours  of  over-all  sound  pressure  level 

have  been  plotted  in  Figs-  23  and  24.  Also  shown  are  the  corresponding 

-55- 


measured  sound  pressure  levels  given  in  Ref.  13*  There  is  good  agreement 
at  MRT  but  the  calculated  values  at  Max  A/B  are  about  5db  higher  than  the 
measured  values.  A  measured  octave  band  frequency  analysis  was  available 
from  Ref.  13  for  a  location  50  ft.  from  the  nozzles  at  ^0°  from  the  center 
line.  This  and  the  corresponding  calculated  values  have  been  plotted  for 
comparison  in  Fig.  25. 

The  5&b  discrepancy  at  Max  A/B  has  not  yet  been  fully  explained. 

Since  Fitch's  method  is  a  modification  of  an  earlier  method  by  Franken 
(Ref.  U2)  a  calculation  was  made  by  Franken fs  method.  Without  afterburner 
the  only  difference  is  that  Franken  in  effect  uses  a  value  of  8  for  n  in 
Item  2.  Since  this  adder  was  negligible  regardless  of  the  value  of  n  the 
two  methods  gave  identical  results  at  MRT.  At  Max  A/B  Franken  uses  the 
same  velocity  correction  as  that  for  MRT,  still  negligible.  The  correc¬ 
tion  for  the  addition  of  the  afterburner  is  a  function  of  the  thrust  of 
the  turbojet  at  100$  rpm.  This  adder  is  obtained  from  a  curve  and  is 
applied  to  all  the  contours.  In  this  case  the  adder  was  8db  instead  of 
the  values  ranging  from  7db  to  12db  for  the  different  contours  in  Item  2. 
The  angular  rotation  is  the  same  as  in  Fitch’s  method.  The  results  are 
plotted  in  Fig.  26.  The  calculation  is  still  approximately  5db  higher 
than  the  measurements  at  most  locations .  However,  concurrent  with  this 
study,  research  into  noise  generation  by  high  temperature  jets  has  been 
carried  out  under  the  Phase  I  program  of  this  contract.  Results  indicate 
that  for  high  temperature  (afterburning)  Jets  estimates  of  the  noise 
generated  cannot  be  immediately  extrapolated  from  lower  temperature  jet 
noise  data  by  means  of  a  simple  velocity  or  thrust  relationship.  It 
appears  that  further  investigation  into  the  generation  of  noise  by  high 
temperature  jets  is  warranted,  and  th.*r.  subsequent  improvement  of  noise 

prediction  techniques  is  required. 
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AXIAL  DISTANCE  FROM  JET  NOZZLE  -  X/D 


30 


i 

i 

I 


Figure  20.  Near  Sound  Field  of  a  Turbojet  (Ref.  39) 
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JET  VELOCITY,  V  IN  FT/SEC 

Figure  21  Velocity  Dependence  of  Over-all  SPL  Contour  Pattern  (r«.£ 
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Figure  22.  Spectrujn  of  Jet  Noise  Near  Field  (Ref.  39) 
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3.  Random  Noise  and  Vibration 


A  random  variable,  such  as  jet  noise,  is  varying  in  time  such  that  its 
instantaneous  values  at  prescribed  times  are  not  predictable  and  there  are 
no  periodicities.  Fig*  27  shows  samples  of  a  single  random  variable  observed 
at  different  time  intervals.  Presumably  one  could  continue  to  sample  this 
function  indefinitely  and  obtain  different  curves.  There  are,  however, 
statistical  parameters  of  such  a  random  function  which  are  constant  under 
certain  frequently  encountered  conditions. 

Let  us  consider  a  specific  interval  x  to  x+Ax  of  the  random  variable  and 
find  the  total  time  that  the  variable  lies  in  this  interval  for  all  the 
available  samples.  We  then  define  the  amplitude  probability 

P  (x,  x+Ax)  «  ZAt^  (6) 

T 


where 

At^  =  time  during  which  x^  is  in  interval  between  x  and  x+Ax 
T  =  total  time  over  which  phenomenon  is  studied 


The  amplitude  probability  density  is 

p(x)  -  L  P(x,  x+Ax)  (7) 

Ax  0 


Assume  that  we  have  a  large  assembly  or  ensemble  of  random  time  records. 
Let  N  be  the  number  of  members  of  the  ensemble.  At  any  given  time  such  as 
t  *=  t  in  Fig.  27  the  x(tp  will  be  different  in  the  various  records  but  we 
can  count  the  fraction  of  records  n/N  for  which  x(t^)  is  equal  to  or  less 
than  a  specified  x.  This  fraction  is  plotted  versus  x  and  tends  toward  a 
fixed  curve  as  N  approaches  oo .  This  curve,  shown  below,  is  called  the 
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probability  distribution  P(x,  t^)  of  the  random  function  at  time  t^. 


If  the  derivative  exists,  we  define  the  probability  density  function  as 

p(x,  tx>  =  -|r  p(x’  ti)  <8) 

so  that  the  probability  of  x(t^)  having  a  value  in  the  interval  between  x 
and  x+dx  is 


dPfx.tp  =  p(x,  t^)dx  ■  P(x+dx,  tj)  -  PCx.tj) 


(9) 


Also 


P(x,tj) 


p(x,t1)dx 


-  oo 


(10) 


A  random  process  is  stationary ; if  any  translation  of  the  time  origin 
Teaves  the  statistical  properties  unaffected.  It  is  ergo die  if  each  sample 
is  representative  of  the  process*  In  particular  if  the  statistical  proper¬ 
ties  are  not  changing  as  a  function  of  time,  the  probability  functions  at 
any  other  time  should  be  the  same  as  that  found  at  t  B  t^.  Also  we  should 
get  the  same  result  by  counting  or  averaging  over  a  time  interval  of 
sufficient  length  as  by  counting  or  averaging  over  the  ensemble  of  functions. 

The  most  commonly  encountered  random  process  has  a  Gaussian  or  "normal" 
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probability  density  distribution. 

.  .ft.-..*?? 
“  ~dtfsr  e  2  °2 


(ii) 


First  and  Second  Moments 

The  first  moment  of  a  random  function  la  its  average  value  defined  by 


-oo 

x  •  /  x  pCx.tjJdx 
'oc 

The  second  moment  Is  defined  by 

_  f~ 

x2  *  /  x2  pfx.tjjdx 
-  oo 


(12) 


(13) 


The  square  of  the  standard  deviation,  or  variance  Is  defined  by 


o2  »  /  (x  -  x)2  p(x,t1)dx  (14) 

Joo  roc 

Note  that  x  Is  a  constant  and  that  /  p(x, t^)  dx  ■  1  since  p(x,t^)  Is  a 

probability  density  function.  * oc 

Then 


2 


a 


2 

x 


(15) 


If  the  process  has  zero  mean  as  is  often  the  case  the  variance  Is  equal  to 
the  mean  square  value. 

For  a  stationary  and  ergodic  random  process  the  moments  are  independent 
of  time  and  time  averages  can  be  used  instead  of  the  ensemble  averages. 

Thus  the  average  value  of  x  is 


<  x  >  =  L 

T-*oo 

» I 

x(t)  dt 

(16) 

<r  x2>-  J 

7  -*>oo 

x2(t)  dt 

(17) 

-67- 


Pover  Spectral  Density 


Random  vibration  may  be  considered  as  the  sum  of  a  large  number  (tending 
to  infinity)  of  harmonic  vibrations  of  appropriate  amplitude  and  phase. 


Consider  a  random  function  in  the  time  interval  2T.  It  can  be  represented 
by  the  real  part  of  the  Fourier  series, 

x(t)  £  C  e  jnV  «  C  +  211.1c  lcos(no)  t  -  9  )  (18) 

n  o  nMi  |  n  i  o  n 

f*  E5  -An  1 


where 


co 


o  T 


(19) 


Cn  =  /  x(t)  €  ’Jna)oT  dr 


(20) 


The  mean  square  value  of  x(t)  is 
T 


<  x  >  *  L 

1^00 


3?  /  x2<C>  *  '  U  W  * 


(21) 


-T 


The  contribution  to  the  mean- square  value  in  any  frequency  interval  is  the 

summation  of  all  the  components  within  that  frequency  band. 

If  now  the  interval  2T  is  doubled,  co^  is  one- ha If  its  former  value  and 

the  number  of  spectrum  lines  in  any  frequency  interval  is  twice  the  former 

2 


number.  The  coefficients  C 


must  be  reduced  to  one-half  their  former 
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values  because  the  mean  square  of  the  entire  spectrum  Is  a  constant  for  an 

assumed  efgodlc  process.  It  follows  that  for  any  given  frequency  interval 

*  12  2 
the  sum  of  the  jCn  j  In  that  interval  Is  a  constant,  £<-x  > .  The  ratio 

2 

of  £&x  >  divided  by  the  frequency  interval  Is  called  the  power  spectral 
density  in  that  frequency  interval 


W(oj) 


&-<x2  ^ 

£a> 


As  T  approaches  >c  ,  co>  approaches  0  and  the  spectrum  becomes  continuous 
and  the  power  spectrum  density  is  given  by 


W(u>) 


L 


A<x2> 


d<-x2> 

do) 


(22) 


From  Eq.  22  we  obtain  the  mean  square  value 

.oo 

W(cd)  dm 


<x2> 


/ 


(23) 


The  dimensions  of  the  power  spectral  density  are  the  square  of  x  per 
radian  per  second. 


Input-Output  Relationship 


If  the  input  to  a  linear  system  is  x(t)  then  the  output  is 
t 

-oo 


y(t) 


( 


x(T)h(t-r)dT  *=  /  x(t-|)h(?)d£ 


-  OO  -'o 

where  h(t)  is  the  response  of  the  system  to  a  unit  impulse  8(t), 


(24) 

Assume 


that  the  input  was  0  for  t><0  and  change  the  lower  limit  to  0  in  both 
integrals.  Let  x(t)  =  then  Eq.  24  gives 
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y(t) 


'GO 


^(t~T>h(r)dT  .  «**  c^hMdr  *  e**A(J«)  (25) 


'OO 


where 


A(ja)) 


l'*"  h(T)  dT 


(26) 


is  defined  as  response  function  of  the  system.  y(t)  of  course  is  the  response 
to  a  harmonic  excitation.  Since  a  sinusoidal  input  is  given  by  either  the 

i 

real  or  imaginary  part  of  the  response  to  a  sinusoidal  input  is  given 

by  the  real  or  imaginary  part  of  y(t).  Since  h(*r)  ia  zero  for  t^O  the 

lower  limit  in  Eq.  26  can  be  changed  to  -oo  without  altering  the  result  and 

A(Jcc)  is  then  the  Fourier  transform  of  the  impulsive  response  function  h(T). 

If  the  response  to  an  Impulse  is  known,  Eq.  26  can  be  used  to  obtain  the 

response  function  A(jcu).  Often,  however,  it  may  be  more  convenient  to 

calculate  the  response  to  a  harmonic  excitation  directly. 

Consider  now  the  input  function  given  by  Eq.  18  which  is  a  random 

function  over  a  finite  time.  Substituting  in  Eq.  24  we  obtain 

roo 

y(t) 


I  OO  JWftj(t-T) 

/  h(T)  2  C  €  dr 

J  n  «  -c«  n 


-co 


,SL 


OO 
2 

n  “  -oo 


C  e 
n 


jncD0t 


h(r)  e 


■jno) 


T  1 


dr 


.(L  f 


n  «  -oo 


C  € 
n 


A(jnu>0) 


Note  that  here  we  have  a  summation  of  terms  like  the  one  given  in  Eq.  25  and 
26  for  a  single  harmonic.  Taking  the  real  part  we  have 


V  OO 

y(t)  =  C  A(o)  +  2  £ 

n  -  1 


|A(jnu)o) 


cos(n^  t  +  a  J 
o  n 


(27) 
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Ftom  Eq.  16  the  mean- square  value  of  the  response  :1s 


L 

T- 


►oo 


1 

2T 


rT 

I  2.. 

y  (t) 


J 

-T 


dt 


Cq2A2(0) 


oo 
+  2Z 
n  ■  1 


A(jncJ)o) 


We  find  the  limit  by  the  same  method  used  in  connection  with  Kq.  21 
power  spectral  density  of  the  response  is 


W  (tO)  «  W  (0>) 

r  s 


A(ju>) 


and  the  mean* square  output  is 

-oo 


<y2^ 


(co)  do) 


.00 

Wg(co) 


A(jcu)  da> 


Response  of  Single  Degree  of  Freedom  System 


The  differential  equation  for  the  steady  state  displacement  of 
in  the  simple  mechanical  oscillator  shown  is 


mjr  +  cy  +  ky*F€ 


jcut 


Let 


k 


m 


5 


C 

hm 


(28) 


The 


(29) 


(30) 


the  mass 


(31) 
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Then  Eq.  31  may  be  written  in  the  dimensionless  form 


(32) 


The  displacement  response  function  for  unit  applied  force  is 


(33) 


and 


Eq*  33  and  34  give  the  response  and  the  square  of  the  response  to  a  unit 
sinusoidal  excitation  force* 

When  the  damping  is  small  the  bandwidth,  26o^,  of  the  response  of  the 
mechanical  oscillator  is  small.  If  the  power  spectral  density  of  the  input 
function  is  essentially  constant  in  the  neighborhood  of  the  resonant 
frequency,  co  the  power  spectral  density  of  the  response  is  approximately 

Wr(co)  -  Wa(o>n)  Ia(Jo)  2  (35) 


where  W  (c&)  has  the  constant  value  W  (cj  ),  its  value  at  cu  *  ca  .  The  curve 
s  s  n  '  n 

of  the  response  power  spectral  density  versus  o>  will  be  very  similar  to  the 
curve  of  the  square  of  response  to  sinusoidal  excitation.  Under  the  same 
conditions,  the  mean  square  response  is 
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This  integral  may  be  evaluated  by  convolution  in  the  complex  plane  to  give 


^  y2^ 


W  (to  ) 
sN  n 


l 


(36) 


The  response  of  a  lightly  damped  (8<0*1)  mechanical  oscillator  to 
random  excitation  has  the  characteristics  of  a  sinusoidal  wave  with  slowly 
varying  amplitude  and  phase*  Miles  (Re£*43)explained  this  qualitatively  by 
noting  that  such  a  system  acts  as  a  narrow  band  filter  passing  only  those 
frequencies  in  the  vicinity  of  Then  by  analogy  with  the  phenomenon  of 

"beats11  the  sum  of  two  harmonic  waves  of  approximately  the  same  frequency  is 
replaced  by  a  wave  having  a  frequency  equal  to  the  average  of  the  two  waves 
and  an  amplitude  envelope  that  fluctuates  at  a  rate  equal  to  the  difference 
frequency.  Thus  it  can  be  inferred  that  the  sum  of  all  the  frequencies  in 
the  pass  band  is  a  sinusoidal  wave  with  frequency  and  an  amplitude  that 

exhibits  a  random  fluctuation  at  a  rate  equal  to  the  bandwidth,  280^, 
or  less. 

It  has  been  stated  by  many  writers  that  the  probability  distribution  of 
the  envelope  is  a  Rayleigh  probability  density  function.  Davenport  and  Root 
(Ref.  4£),for  example,  write  the  response  function  of  a  narrow-band  Gaussian 
process  as 

y(t)  =  V(t)  cos  [oi  t  (37) 

in  j 
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where  th$  variations  of  V(t)  and  <5(t)  are  slow  compared  to  those  of  cos  om:* 
The  function  cart  be  expressed  as  a  Fourier  series  over  the  interval 


oo 

y(t)  ■  £  (y  cos  no>  t  +  y  sin  ncu  t) 

-  w  cn  o  J  sn  o 

n  »  1 


where  cd  ® 
o 


2% 

T 


cn 


r  T 
2  / 

-y  I  y(t)  cos  na}Qt  at 


sn 


-L 

T 


y(t)  sin  nco  t  dt 
o 


(38) 


(39a) 

(39b) 


These  coefficients  are  Gaussian  random  variables  which  become  uncorrelated 
as  T~>oo , 

The  mean  frequency  of  the  resonant  response  is  introduced  by  writing  no)Q  in 

Eq.  38  as  (na>  -  a>  )  +  o  ,  where  a>  =  2k£  ,  and  expanding  the  sine  and  cosine 
o  n  n  n  n 

factors.  The  result  is 


y(t)  ■  y.(t)  cos  a  t  -  y  (t)  sin  ffl  t 
c  n  o  n 


(40) 


with 


00 

yc(t)  -  2 

C  n  -  1 


ycn  C08  <n%  *  an>t  +•  ysn  8ln  (rWo  '  “n)C 


(41b) 


,  ,  00 

VAt)  -  £ 

8  n  -  1 


ycn  8in  (n“o  ‘  wn)t:  '  ysn  C08  (nQo  ' 


(41b) 


Then,  expanding  Eq.  3  7  and  comparing  with  Eq.  4(^  we  see  that 


yc(t)  -  V(t)  cos  4(t) 


yfi(t)  ■  V(t)  sin  ^(t) 


(42a) 

(«b) 
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It  follows  that 


V(t) 


V 


y„2oo 


+  y.  <t) 


(43a) 


Ht) 


tan 


-1 


ys(t) 

yc(t) 


(43b) 


The  only  non- vanishing  terms  in  the  sums  of  Eq.  41  are  the  ones  for  which 

the  values  of  nu>  fall  in  the  response  band  near  o>  .  Therefore  y  (t)  and 
o  n  c 

y  (t)  have  frequency  components  only  in  a  narrow  ban*-  centered  on  zero 
6 

frequency.  The  envelope  and  phase  likewise  have  components  in  the  neighbor¬ 
hood  of  zero  frequency. 

Let  y  and  y  refer  to  the  possible  values  of  y  (t)  and  y  (t) 

Ct  St  C  8 

respectively.  They  are  the  sums  of  independent  Gaussian  random  variables 
having  zero  means 


E(yct)  *  E(yat)  »  0  (44) 

Here  E  (  )  stands  for  the  expected  value  in  a  statistical  sense.  The  mean 
square  of  y  is  obtained  from  Eq.  41a. 


e<4> 


oo  o* 

-  Z  Z 
n»l  m*! 


+ 


+ 


+ 


E(ycn  W  cos  (nV“n)t:  cos  (niwo  '  "nJt 

E(ycn  ysm)  cos  (nV®n>l:  8ln  (lD"o  “  ^ 

E(ysn  ycm}  Bin  cos  (mo>0  -  0R)t 

E(ysn  ysm)  8in  <na>0'Vt  8in  (®"o  *  V* 


-75- 


Since  the  coefficients  are  independent 


E(y  „ 

y  ) 

**  0 

for 

cn 

cm 

E(ye 

y  ) 

=  0 

for  m  ^  n 

sn 

sm 

E(y 

wcn 

=  0 

for  all  m  and  n 

E(y 

7cn 

^CID^ 

■  E(y  y  ) 
sn  sm 

2 

«  E(y  )  for  m  ° 
wcn 

and 


2  oo  2 

E  (y  )  =  £  E(y  ) 

wct  ,  ycn 

n  «  1 


2  2 
cos  (nto^-to) t  +  sin  (nci^-ai) t 


E(y^) 


(45) 


Similarly  it  may  be  shown  that 


E<V>  “ 


(4  6) 


The  variance  of  the  coefficients  is 


a2(yct)  =  E(y2t)  -  E(yct> 


02(yflt)  «  E(y2t)  -  E(yst) 


Substituting  from  Eqs*  44,  45,  and  46,  we  obtain 


a2(yct)  -  s2(yet)  =  E(y2) 


Let  cr  =  o(y).  The  covariance  of  y  and  y  ^  is  obtained  from  Eqs.  41a 
y  Jt  ct  'st 

and  ^b. 
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E<?ct 


n«l  m*l 


E(y 

-  E(y 
+  E(y 

-  E(y 


cn 


cn 


sn 


an 


y  )  cos  (net  -03  )t  sin 
J  cm  on 


y  )  cos  (nco  -cd  )t  cos 
J  am'  v  o  n' 


y  )  sin  (na>  -to  )t  sin 
^car  v  o  ny 


y  }  sin  (no)  -cd  )t  cos 
J  set  v  o  n 


<nV“n}t 

(ncDo-con)t 

<WDo"®n>t 

(nc“o'“n)t 


This  reduces  to 


E<yct  yet) 


oc  2 

S.  B<S'cn> 

n*=l 


cos(n<Do-cDn)t  8in(ncoo-a^)t 

-sin(na)  -cd  )t  cos(na)  -co  )t 
o  n  o  n 


The  random  variables  y 

ct 

zero  means  and  variances 

y 


and  y  have  been  shown  to  be  independent  with 
st 

Their  probability  density  then  is 


P(yct) 


€Xp 


P^st*  "  sf&o  ***  l~8t 


and  their  joint  probability  density  is 


P(yct» 

yst  )  =  P(yct)  •  P(yst) 


2n  o 


exp  - 


2  2 
u  +  y 

'  c t  yst 


2a 


(47) 


The  joint  probability  density  function  of  the  envelope  and  phase 
variablea  and  6^  can  be  found  from  Eq.  47  by  a  simple  transformation  of 
variables. 
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p(Vt,  it)  -  p(yct,  y8t)  |j| 

The  Jacobian  of  the  transformation  equations,  Eqs*  42a  and  42b,  is  given  by 


COS  4>t 

sin 

S' 

CO 

rt 

-V£  sin 

V  cos  i t 

Then  with  from  Eq.  43a 


p(vt>  ^t) 


(48) 


Since  Eqs.  42a  and  42b  define  a  transformation  only  for  2r  0  and 

0-  ^  2jt  the  joint  probability  of  Eq.  48  applies  only  over  these  intervals 

and  is  zero  elsewhere. 

The  probability  density  function  of  the  envelope  alone  is  obtained  by 
integrating  Eq,  48  with  respect  to  over  the  range  0  to  2rt, 


P(Vt) 


for  V  ~  0 


(49) 


This  is  the  Rayleigh  probability  density  function. 
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FatigUe  of  metals  is  caused  by  repeated  plastic  deformation  during 
cyclic  loading*  Cracks  are  initiated  at  stress  concentrations  which  may 
be  either  geometric,  for  example,  at  points  of  attaching  a  panel  to  a  rib 
or  near  a  hole,  or  they  may  be  dynamic,  a  function  of  the  mode  of  vibration* 
The  stress  life  curve  developed  for  constant  amplitude  cycling  of  the  load 
has  for  many  years  provided  the  basic  data  for  designing  structures  to 
prevent  fatigue  failure.  This  curve  is  referred  to  as  the  S-N  curve,  S 
for  stress  and  N  for  cycles  of  life.  In  acoustic  fatigue  a  problem 
arises  because  of  the  high  rate  at  which  cycles  of  load  are  applied  rather 
than  a  high  magnitude  of  stress.  Even  so  the  design  of  a  given  structure 
may  be  limited  by  the  acoustic  load  and  it  is  desirable  to  use  the  highest 
possible  stress  to  increase  the  payload. 

The  pressure  loads  produced  by  jet  noise  or  boundary  layer  noise  are 
random  functions  of  time.  The  amplitude  distributions  are  approximated 
quite  well  by  the  “normal"  or  Gaussian  distribution  function.  With  this 
distribution  the  major  part  of  the  damage  to  a  structure  is  caused  by  a 
relatively  few  cycles  at  the  higher  stress  levels.  Collection  of  fatigue 
life  data  with  random  excitation  is  therefore  relatively  inefficient  and 
time  consuming.  Although  some  life  data  are  available  for  random  excitation, 
a  much  greater  amount  has  been  taken  for  S-N  curves  at  constant  amplitude 
sinusoidal  excitation.  It  is  desirable  therefore  to  utilize  this  vast 
experience  if  possible  by  finding  a  relationship  between  cumulative 
fatigue  damage  under  random  loading  and  available  S-N  data. 

Much  work  is  being  done  on  life  testing  of  both  samples  and  structures 
utilizing  random  noise  sources  as  well  as  programmed  sinusoidal  testing 

to  duplicate  the  essential  effects  of  random  noise. 
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4A.  Linear  Accumulation  of  Damage 

Probably  the  roost  widely  used  method  for  analyzing  fatigue  damage  and 
predicting  structural  life  is  Miner's  Rule.  It  is  relatively  simple  to 
apply  and  gives  reasonable  results.  Numerous  other  methods  have  been 
proposed  which  are  aimed  at  improving  the  accuracy  by  including  the  effects 
of  stress  interaction  and  non-linearities.  Some  of  these  methods  will  be 
discussed  in  later  paragraphs. 

Miner's  method  has  the  advantage  that  it  makes  use  of  the  available 
S-N  data  and  provides  for  the  linear  accumulation  of  damage  without  regard 
to  time  history  of  loading.  The  order  in  which  variable  amplitudes  of 
stress  or  rest  periods  are  applied  are  not  included.  It  is  assumed  that 
the  stress  cycles  are  fully  reversed.  This  is  approximately  true  for  a 
lightly  damped  structure  vibrating  in  a  single  predominant  mode  with  zero 
mean  stress  even  though  the  excitation  is  Gaussian. 

Miner's  Rule  is  based  on  an  S-N  curve  such  as  Fig.  28.  Let  the  load 
be  cycled. for  n^  cycles  at  stress  amplitude  9^. 

Stress 


Fig.  28 


Since  the  total  life  at  this  stress  is  N.  cycles  it  is  postulated  that 

^1 

there  has  been  damage  equivalent  to  an  rr*  proportion  of  the  expected  life. 

1 

It  is  clear  that  failure  should  occur  when  this  proportion  equals  1.  If 
the  load  is  also  cycled  at  S2  for  n*  cycles  the  total  damage  is 

nl  +  n2 


N. 


N, 
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Thus  for  any  number  of  stress  levels  the  condition  for  failure  is 


+ 


+  ^i 

N, 


1 


Let  the  total  number  of  cycles  at  the  various  loads  be  given  by 

N  **  Sni 


Then  the  condition  for  failure  can  be  written 

nl  +  n2  +  ...  +  ni  -  1 

NHj  N  K2  N  N£  N 


For  a  random  load  the  probability  of  occurrence  of  a  given  stress 
amplitude  is 


Hence 


pi 


+  pi 


i 

N 


or  the  total  life,  in  terms  of  cycles  to  failure,  is  given  by 
N  =  1 


(51) 


For  a  Gaussian  distribution  of  stress  the  distribution  of  fhe  peak 
amplitudes  in  the  response  of  a  single  degree  of  freedom  vibrating  system 
is  given  by  the  well-known  Rayleigh  distribution  function.  (See  Eq.  49) 


p(a)  =  “a  e 


-  (  s2/2 a2 


(52) 


v*iere  s  is  the  magnitude  of  peak  stress  and  a  is  the  root  mean  square 
stress  for  the  case  of  zero  average  stress.  Since  this  is  a  continuous 
function  of  s  the  summation  in  Eq.  51  can  be  replaced  by  an  integral. 
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Then  the  number  of  random  cycles  is 


where  N(s)  is  expressed  by  the  equation  for  the  constant  amplitude  S-N 
curve. 

Miles  (Ref.  43)  used  the  function 

"Hrf 

where  s^  is  the  stress  for  which  the  equation  would  predict  failure  at 
one  cycle  although  the  equation  should  not  be  extrapolated  this  far.  a  is 
the  reciprocal  of  the  slope  of  a  log  S  vs.  log  N  curve.  Miles  derived  an 
expression  for  the  equivalent  stress,  or  "reduced  stress",  which  would 
produce  the  same  fatigue  damage  as  the  random  load  after  the  same  total 
number  of  cycles  of  loading. 


Shanley  (Ref.  45)  has  suggested  that  a  should  be  replaced  by  2a  in  this 
equation.  Belcher,  et  al  (Ref.  40)  have  shown  that  the  log  S  -  log  K 
curve  may  lead  to  highly  conservative  results.  They  have  instead  obtained 
graphical  or  numerical  solutions  of  Eq.  53  using  data  from  actual  S-N 
curves  for  various  materials.  The  procedure  as  explained  below  is  based 
on  an  example  including  a  table  and  curves  from  Ref.  39. 


From  Eqs.  52  and  53  we  obtain 
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or  for  numerical  calculation 


1  e 


_I  (s/c)' 
-2 


ii) 


N(s) 


I 


>(i)  4) 

N(s) 


(54) 


Having  an  S-N  curve  for  the  material,  Fig.  29, and  the  Rayleigh  probability 
distribution  curve  Fig.  30,  we  can  solve  Eq.  54  for  N^  for  any  selected 
value  of  the  random  rms  stress, c  •  By  repeating  the  calculations  for  other 
values  of  a  we  obtain  the  points  for  plotting  a  random  S-N  curve. 

The  calculation  for  each  o  is  facilitated  by  a  table  such  as  Table 
XVII.  First,  a  value  for  the  rms  stress,  c  ,  is  assumed,  then  values  of  the 
relative  stress  —  are  chosen  at  discrete  intervals.  For  each  relative 


S  I « 

stress,—  ,  the  relative  number  of  cycles, 


is  read  from  the  Rayleigh  dis¬ 


tribution  curve,  Fig.  30.  Multiplying  each  relative  stress  by  o  gives  the 

corresponding  stress,  s,  which  can  be  used  to  enter  the  S-N  curve,  F.ig.  29, 

to  obtain  N(s),  the  number  of  allowable  cycles  at  that  stress.  The  ratio 

P(a) 

— is  the  relative  damage  or  damage  density  and  is  plotted  versus  ~ 

N  (8)  0 

in  Fig*  31*  At  the  maximum  damage  density  we  can  read  the  relative  stress 

and  obtain  the  peak  damage  stress,  s  ,  *  The  area  under  this  curve  repre- 

pd 

sents  the  summation  indicated*  in  Eqa  54  and  can  be  computed  by  multiplying 

i&i 


the  sum  of 


N(s) 


as  givan  in  the  last  column  in  Table  XVII  by  the  relative 


stress  interval  •  The  result  Is  the  reciprocal  of  the  predicted  life 

cycles  at  the  assumed  rms  level  of  the  random  stress*  The  same  process 
was  repeated  for  other  vabes  of  o  and  the  results  plotted  in  Fig*  29  as  a 
random  S-N  curve.  The  peak  damage  stress,  Spj ,  is  also  plotted.  It  is 
seen  that  s  ^  is  3  to  4  times  the  rms  stress*  Referring  to  the  Rayleigh 
distribution  curve  we  see  that  it  is  only  one  per  cent  or  so  of  the  stress 

i 

cycles  which  cause  most  of  the  damage.  In  a  siren  test  the  stress  is 
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igure  29.  Comparison  Between  Actual  and  Predicted  Spectrum  Life  (Ref>  39) 


RAYLEIGH  PROBABILITY  CURVE 
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CUMULATIVE  DAMAGE  TABLE  FOR  ANNEALED  TITANIUM 
0=15, 000  PSI  (ASSUMED) 
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^  P(oj  =  0.2  x  27.4  x  10"6  =  5.49  x  10"6 
^  N(s)  ! 

N  =  — - — -  =  182,000 

5.49  x  10* 6 

TABLE  17.  CUMULATIVE  DAMAGE  TABLE 
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Figure  ii.  Damage  Density  Curve  (Ref  39) 
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Thus  the  siren  can  apply  enough  cycles 


usually  pretty  close  to  8p{j* 
to  fail  the  sample  in  a  few  minutes  whereas  the  equivalent  random  source 
of  noise  would  require  hours. 

Random  S-N  curves  for  several  materials  as  well  as  honeycomb  struct¬ 
ures  have  been  computed  and  are  given  in  Ref.  39. 

There  has  been  criticism  that  Minerls  rule  ignores  factors  which  are 
known  to  influence  fatigue  life.  For  example,  the  time  sequence  of  variable 
loads  is  neglected.  In  two  level  block  testing  on  mild  steel  samples 
Kommers  (Ref.  46)  found  that  improved  fatigue  strength  could  be  obtained 
by  applying  the  low  loads  first.  The  opposite  effect  was  found  in 
aluminum,  (Ref.  47).  When  many  blocks  of  stresses  alternating  between  low 
and  high  are  used  the  sequence  effect  becomes  less  important  as  the  number 
of  blocks  is  Increased  but  it  apparently  does  not  completely  disappear. 

Miner's  rule  makes  no  allowance  for  rest  periods  and  there  is  no  way 
to  include  the  effects  of  static  loading.  Residual  tensile  stress  is 
detrimental  but  compressive  stress  may  be  beneficial, 

Schjelderup  (Ref.  48)  has  analyzed  the  strain  gage  response  of  sheet 
rib  structures  to  excitation  by  a  jet  engine  on  an  aircraft  and  in  a  test 
cell.  He  found  that  the  Rayleigh  distribution  accurately  described  the 
distribution  of  peaks  in  both  cases.  Thus  for  a  lightly  damped  single  degree  v. 
of  freedom  system  the  response  is  a  sinusoid  with  varying  amplitude.  But 
with  many  modes  the  response  is  a  complex  wave  form.  According  to 
Schjelderup  some  distribution  other  than  the  Rayleigh  must  be  used  for 
multimode  response  but  he  did  not  offer  a  substitute.  As  we  shall  see 
later  Belcher  (Ref.  40)  applies  a  multimode  correction  factor  in  his  method. 

Many  writers  have  reported  that\"*^i  can  often  be  larger  or  smaller 
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than  unity  at  failure.  For  example  Fralich  (Ref,  49)  has  reported  that 
the  fatigue  life  was  over  estimated  for  some  notched  SAE  4130  steel  beam 
samples.  Ref.  47  cites  the  report  by  Hooke  and  Head  that  Miner's  rule 
over  estimated  the  life  of  notched  samples  by  as  much  as  3-1/2  times  at 
low  stresses  and  somewhat  less  at  higher  stresses.  The  same  reference 
also  cites  the  work  of  R,  W*  Fralich  on  notched  test  pieces  of  7073  AL 
alloy  which  were  subjected  to  reverse  bending  within  the  range  of  their 
natural  frequency.  Contrary  to  the  results  of  Hooke  and  Head  he  found  that 
Miner's  rule  under  estimated  the  life  at  high  stresses  and  predicted  it 
accurately  at  low  stresses.  Ref.  47  mentioned  that  a  program  of  research 
was  undertaken  at  the  University  of  Southampton  to  investigate  the  con¬ 
tradictions  in  the  results  reported  in  random  loading  studies.  A  preliminary 
report  by  M.  T.  Lowcock  and  T.  R.  G.  Williams  concluded  that  Ml)  an 
accurate  estimate  of  the  fatigue  life  under  variable  loading  conditions 
can  be  obtained  from  the  rms  of  the  stress  and  the  S/N  curve  of  the 

material.  2)  There  is  a  discontinuity  in  the  S/N  curve  of  L.73  alloy  at 

2 

a  stress  level  of  16  tons/in.  ,  and  the  shape  of  the  variable  loading 
fatigue  graph  is  modified  by  the  discontinuity." 

4B.  Freudenthal  Method 

The  discontinuity  in  the  S-N  curve  has  been  investigated  by  Freudenthal. 
(Ref.  50)»  Based  cn  tests  on  small  plain  samples  there  appear  to  be  two 
distinct  mechanisms  of  fatigue.  The  first  is  called  high  level  fatigue. 

It  occurs  at  either  high  stress  or  low  frequencies  such  that  the  deformation 
I 8  akin  to  that  in  unidirectional  "static"  deformation.  Failure  is  related 
to  a  critical  amount  of  total  plastic  deformation  or  of  strain  hardening 
due  to  distortion  and  breakdown  of  the  grain  structure.  It  corresponds  to 
the  short-life  portion  of  the  S-N  curve  which  has  a  high  slope* 
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Rapid  cycling  stresses  of  small  amplitudes  produce  less  strain  harden¬ 
ing  and  no  significant  distortion.  Instead  a  multitude  of  fine  slip  bands 
are  concentrated  in  striations  so  that  individual  dislocations  cannot  be 
resolved.  This  corresponds  to  the  longer- life  portion  of  the  S-N  curve 
with  smaller  slope. 

The  two  mechanisms  differ  in  the  way  energy  is  released  as  heat  as 
well  as  the  way  cracks  are  initiated.  The  transition  between  the  two 
mechanisms  should  be  sharp  but  its  location  depends  on  the  frequency  of 
cyclic  stressing,  the  initial  condition  of  the  grain  structure,  the  blocking 
of  the  primary  slip  system  by  dislocations  and  the  intensity  of  the  mean 
stress. 

With  varying  stress  amplitudes  there  is  an  interaction  between  the 
two  mechanisms.  In  particular  the  high  level  mechanism  has  the  effect  of 
accelerating  damage  by  the  low  level  mechanism.  Thus  a  relatively  few 
cycles  at  the  high  level  will  shorten  the  fatigue  life  under  the  low  level 
stress  amplitudes  far  beyond  any  damage  immediately  associated  with  the 
high  level  cycles*  Freudenthal  considered  using  an  interaction  factor  to 
predict  this  effect.  He  concluded  however  that  it  would  be  better  to 
create  a  ficticious  S-N  curve  for  variable  loading.  He  defined  the  S-N 
curve  in  two  segments  as  follows; 


Nt<  N 
N  >  N 


where  (N,  S)  is  the  point  on  the  S-N  curve  at  the  boundary  between  high 
level  and  low  level  fatigue  and  (N^,  S^)  is  any  other  point  on  the 
ficticious  S-N  curve.  The  stress  interaction  is  introduced  by  changing  the 
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exponent  from  v  to  p.  For  materials  such  as  2024  and  7075  aluminum  and 
4340  steel  8<v<16  and  4<p<8. 

Freudenthal  states  that  the  effect  may  be  different  for  larger  or 
more  complex  samples.  Thus  a  certain  amount  of  testing  would  be  needed  to 
establish  the  parameters. 

4C.  Smith's  Cumulative  Damage  Method 

Ref.  39  outlines  a  method  proposed  by  C.  R.  Smith  which  would  modify 
the  S-N  curve  so  that  Miner's  method  could  be  used*  He  argues  that  the 
principal  errors  in  fatigue  life  prediction  are  due  to  the  residual  stresses 
at  concentrations.  Thus  the  few  cycles  of  high  stresses  in  random  loading 
alter  the  S-N  curve.  The  S-N  curve  can  be  corrected  by  applying  a  preload 
that  is  equivalent  to  the  highest  probable  load  that  can  be  expected  in 
the  first  ten  per  cent  of  service  life.  Although  Miner's  method  can  be 
used  with  the  modified  S-N  curves  new  data  must  be  acquired  for  each  material 
and  load  spectrum*  The  available  S-N  data  cannot  be  used  with  this  method. 

4D.  <3  Methods  of  Fatigue  Analysis 

Fuller  (Ref.  51)  has  developed  a  method  for  predicting  the  life  of  a 
structure  exposed  to  narrow  band  random  stress  as  well  as  broad  band  random 
stress.  It  is  also  applicable  to  variable-cycle  loading . provided  the  stress 
patterns  or  blocks  of  stresses  are  repeated  enough  times  to  make  the  order 
of  loading  unimportant. 

The  order  in  which  stresses  are  applied  affects  the  total  cycles  to 
failure.  For  two  levels  of  stress  cycles  with  greater  then  the 
accumulated  cycle  ratio^^  is  less  than  1.0  at  failure,  while  with  S^<; 
the  cycle  ratio  is  greater  than  1.0.  If  an  S-N  curve  is  determined  using 
samples  which  have  all  been  cycled  at  a  constant  maximum  stress  for  an 
equal  number  of  cycles,  the  new  S-N  curve  will  cross  the  original  S-N  curve 
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at  the  value  of  the  pre-strens.  The  number  of  cycles  at  pre-stress  are  included 
in  the  new  plot. 

Samples  stressed  below  the  pre-stress  will  fail  earlier  than  predicted  by  the 
original  S-N  curve  because  the  pre-stress  cycles  are  more  damaging  than  an  equal 
number  of  cycles  at  the  test  stress.  Furthermore  Freudenthal  has  shown  that  there 
is  a  stress  interaction  which  may  cause  more  rapid  accumulation  of  damage  at  the 
test  stress  after  a  higher  pre-stresB.  If  the  test  stress  is  larger  than  the  pre¬ 
stress  the  total  life  is  generally  longer  than  would  be  indicated  by  the  original 
S-N  curve  at  the  test  stress*  It  Is  thought  that  understr easing  tends  to  relieve 
the  peak  residual  stresses  in  the  sample  and  increase  the  life.  The  new  S-N  curve 
for  the  pre-stressed  namples  would  thus  tend  to  rotate  about  the  pre-stress  point 
on  the  original  S-N  curve.  The  amount  of  rotation  is  approximately  proportional 
to  the  accumulated  cycle  ratio  of  the  pre-stress.  For  no  pre-stress  R=0  and 

there  is  no  rotation  and  all  points  lie  on  the  original  S-N  curve.  For  Rrl.O  the 
sample  must  fail  before  any  test  cycles  are  applied. 

Fig.  32  illustrates  the  rotation  of  the  S-N  curve.  First  s^  is  applied  for 
n^  cycles  and  curve  1-1  is  rotated  Into  2-2.  Then  s2  Is  applied  for  n2  cycles. 

The  curve  is  further  rotated  about  the  point  at  the  intersection  of  the  curve  2-2 
with  the  stress  level  This  results  In  curve  3-3.  It  can  be  seen  that  the 

fatigue  life  will  not  be  less  than  that  given  by  the  original  S-N  curve  for  the 
maximum  stress,  nor  will  it  be  longer  than  that  given  by  the  original  S-N  curve  for 
the  minimum  stress.  The  actual  fatigue  life  for  variable  cycle  loading  would  fall 
in  this  interval  with  its  actual  location  a  function  of  the  load  distribution.  As 
a  practical  matter  the  upper  limit  of  the  interval  is  established  by  extrapolating 
the  highest  slope  tangent  to  the  original  S-N  curve  down  to  the  minimum  stress. 

A  computed  factor,  locates  the  predicted  fatigue  life  as  the  fraction  of 
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figure  32.  EFFECT  OF  PRESTRESS  CYCLES  ON  TOTAL 
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the  interval  from  the  upper  life  limit.  Thus  if  *3  is  near  zero  most  of 
the  stress  cycles  were  at  the  lowest  stress  and  the  fatigue  life  would  be 
at  the  upper  limit  of  the  interval.  If  3=1.0  all  of  the  stress  cycles 
would  be  at  the  maximum  stress  and  failure  would  occur  at  the  lower  limit 
of  the  life  interval  or  at  the  life  shown  on  the  original  S-N  curve. 

For  variable  stress  cycles  With  a  continuous  distribution 
between  a  minimum,  S  ,  and  maximum  ,  Fuller  gives  the  following 
equation  for  3: 


_ 1 _  f  i 

q(SA  '  Sa| 


log  P(S-Sa)  +  q 


ds 


If  the  lowest  stress  is  zero  this  reduces  to 


q  SA 


log  P(S)  +  q 


ds 


q  is  a  stress  sensitivity  factor  which  must  be  selected  by  the  user  of 
this  method  based  either  on  experience  or  on  two- level  block  testing. 

PCS*Sa)  is  the  cumulative  probabilty  distribution  of  the  stress  loading. 

For  a  Rayleigh  distribution  of  peaks  about  zero  mean,  (&  is  simply  a  constant. 


2 

3' 


Fuller  has  analyzed  Trotter's  narrow  band  random  data  (Ref.  52)  using 
the  ^method  and  obtained  very  good  results,  see  Fig.  33.  Q  Is  equal  to 
A/I.  He  also  analyzed  the  results  of  extensive  tests  on  samples  loaded 
with  three  different  broad  band  random  stress  spectra.  With  suitable  choice 
of  the  stress  sensitivity  factor,  q,  the  results  looked  very  promising. 


4E.  Strain  Life  Method 

In  their  analysis  of  sonic  fatigue  life  determined  by  siren  testing, 
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(a)  BETA -METHOD  ANALYSIS  USING 
MAXIMUM  SLOPE  TANGENT 


(b)  BETA- METHOD  ANALYSIS  USING 
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figure  33  COMPARISON  OF  ANALYSIS  AND  TEST  RESULTS, 
TROTTEffe  NARROW  BAND  RANDOM  TESTS  (Rsf.  5i> 
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P.  W.  Smith  and  C.  I,  Malrae  (Ref.  53)  have  used  Miner's  rule  to  determine 
the  fatigue  life  under  random  excitation  from  siren  tests  in  much  the  same 
manner  as  Belcher  (Ref,  40),  They  argue  however  that  better  results  are  ob¬ 
tained  if  strain- life  curves  from  constant  amplitude  testing  are  used  instead 
of  the  more  widely  used  S-N  curves.  There  is  no  difference  of  course  at  low 
stresses  below  the  yield  point  where  strain  is  proportional  to  stress.  Above 
the  yield  point  the  peak  stress  and  peak  strain  are  nonlinearly  related.  The 
question  arises  whether  either  the  local  stress  or  local  strain  can  be  con¬ 
sidered  linearly  related  to  the  overall  structural  response.  The  peak  strains 
usually  are  concentrated  in  localized  areas  such  as  in  the  vicinity  of  a  hole 
or  weld  or  other  point  of  attachment*  The  fatigue  damage  thus  will  occur 
mostly  in  these  small  localized  regions.  It  is  assumed  that  plastic  yielding 
in  these  small  areas  will  not  affect  the  surrounding  material  and  hence  will 
not  alter  the  response  of  the  structure*  Therefore  it  is  assumed  that  the 
local  strains  in  the  fatiguing  area  are  linearly  related  to  the  structural 
response  * 

The  use  of  strain  instead  of  stress  is  important  only  at  short  life¬ 
times  when  a  significant  portion  of  the  damage  is  due  to  strains  above  the 
local  yield  point.  A  constant  stress  lifetime  curve  tends  to  a  lower  slope 
in  this  region  while  the  corresponding  strain-life  curve  continues  to  rise 
smoothly  toward  the  strain  which  will  cause  failure  in  a  ‘'static11  tensile  test. 

Below  the  yield  point  the  constant  stress  S-N  curve  is  easily  converted 
to  a  constant  strain  lifetime  curve.  The  ultimate  strain  for  failure  in  one 
cycle  in  a  static  test  gives  the  extreme  end  of  the  curve.  Until  test  data 
is  accumulated  the  points  in  between  must  be  estimated  by  extrapolation. 

For  high  rms  stresses  in  a  random  test  we  find  that  the  Rayleigh 
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distribution  may  predict  a  small  probability  of  occurrence  of  large  peak 

stresses  exceeding  the  ultimate  strength  of  the  material*  Presumably 

failure  should  occur  by  fracture  of  the  material  much  sooner  than  actually 

occurs  by  fatigue  during  the  test*  Some  workers  have  avoided  this  dilema 

by  arbitrarily  truncating  the  stress  distribution.  An  example  is  given 

in  Ref.  53  which  shows  that  the  use  of  strain  instead  of  stress  gives  a 

more  reasonable  result.  For  a  structure  of  2024  aluminum  without  stress 

-3 

concentrations  an  rms  strain  of  1.9(10)  corresponds  to  rms  stress  of 
20,000  psi.  A  lifetime  of  1.6(10)^  cycles  is  computed  if  the  rare,  very 
high  strains  are  ignored.  The  ultimate  strength  of  the  material^ 64, 000 
psiv  is  only  3.2  times  the  rms  stress.  The  Rayleigh  distribution  predicts 
that  this  stress  will  be  exceeded  about  6  times  in  every  1000  cycles. 
Consequently  the  fracture  of  the  material  would  be  expected  to  be  more 
likely  than  the  accumulation  of  significant  fatigue  damage.  If  strain  is 
used  rather  than  stress  the  problem  of  failure  by  rare  exceedances  dis¬ 
appears.  The  ultimate  strain  for  2024  aluminum  is  about  0.19.  This  is 

-3 

so  much  greater  than  the  rms  strain  of  1.9(10)  that  the  Rayleigh  distri- 

2171 

bution  predicts  it  will  be  exceeded  only  one  in  10  ,  or  practically  not 

at  all. 

4F.  Scatter  of  Fatigue  Data 

It  is  evident  that  there  is  not  yet  a  clear  understanding  of  fatigue 

and  that  there  are  deficiencies  in  the  methods  for  predicting  fatigue  life 

under  variable  loading  conditions.  The  prediction  of  fatigue  life  is 

largely  empirical  and  good  judgement  and  experience  are  essential  to  proper 

evaluation  of  all  the  many  variables  which  enter. 

The  scatter  in  fatigue  data  cannot  be  avoided.  The  formation  of  a 

fatigue  crack  depends  on  residual  stresses  in  the  material  which  arise 
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from  the  forming,  heat  treating,  machining  or  even  handling  of  the  material* 
The  metallurgical  structure  including  the  shape,  size  and  orientation  of 
the  grains  will  affect  the  fatigue*  Variations  are  also  caused  by  a  lack 
of  randomness  in  the  loading,  unknown  stress  concentrations,  stress  inter¬ 
actions,  time  history  of  the  load,  nonlinearity,  pre-loads  as  well  as 
errors  in  measurement.  It  follows  from  the  shape  of  the  S-N  curve  that 
a  change  of  about  5  per  cent  in  the  stress  will  lead  to  about  100  per  cent 
change  in  the  fatigue  life  in  the  long  life  region.  A  somewhat  more 
comforting  view  of  the  situation  is  to  ask  what  design  margin  in  stress  or 
acoustic  load  is  required  to  meet  a  given  design  life.  Eventually  it  may 
be  possible  to  apply  a  statistical  confidence  interval  to  S-N  data.  Fig. 

54  from  Ref.  39  shows  how  the  upper  and  lower  envelopes  of  the  S-N  data 
for  a  much  tested  material  convert  to  random  S-N  curves.  The  conversion 
follows  the  Rayleigh  distribution  of  peaks  and  Miner's  rule  for  damage 
accumulation* 

For  the  development  of  engineering  design  data  it  appears  quite 
reasonable  to  use  Miner1 s  rule  to  take  advantage  of  its  simplicity  and  the 
possibility  of  using  the  large  amounts  of  available  S-N  data.  With 
appropriate  design  margins  it  can  give  good  results.  (Ref.  53  and 
54). 

5.  Stress  Response  to  Acoustic  Loads 

The  acoustic  pressures  in  the  near  field  of  a  jet  engine  produce 
random  vibration  and  stresses  in  the  aircraft  structure.  In  particular 
the  exposed  cover  panels  are  subjected  to  intense  pressures  which  may 
cause  fatigue  failure  although  the  effect  on  stiffening  ribs  or  frames 
may  also  be  important.  Miles  (Ref,  43)  analyzed  the  response  of  a  single 
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degree  of  freedom  system  to  a  random  excitation*  Powell  developed  a  more 
powerful  approach  which  included  several  modes  of  vibration  and  the  space 
correlation  of  the  applied  pressure  field.  (Ref.  55),  Unfortunately  this 
requires  more  information  than  is  usually  available  in  the  design  stage. 

A  knowledge  of  the  mechanical  impedance  of  the  structure  is  needed  which 
depends  on  the  mode  shapes,  the  natural  frequencies  and  the  structural 
and  acoustic  damping  of  each  normal  mode.  Also  the  power  spectral  density 
and  the  narrow  band  spatial  correlation  of  the  pressure  field  are  required. 


5A.  Sine-Random  Stress  Equivalence 

Belcher,  Van  Dyck  and  Eshleman  (Ref.  40)  extended  Mile’s  work  to 

develop  the  equivalence  relation  between  random  and  sinusoidal  stress 

response  so  that  siren  tests  could  be  related  to  random  noise  response. 

The  same  method  is  used  by  Fitch  (Ref,  39).  The  equivalence  relation  can 

be  derived  as  shown  below.  Eq.  34  gives  the  square  of  the  displacement 

response  of  a  single  degree  of  freedom  system  to  a  unit  applied  force. 

If  the  static  displacement  per  unit  force  yo/Fo  is  replaced  by  the  static 

stress  per  unit  pressure,  o,  and  the  result  multiplied  by  the  square  of 

2 

the  rms  pressure,  p  ,  Eq.  34  gives  the  square  of  the  rms  stress  response, 
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Here  the  subscript,  s,  denotes  sinusoidal  test  conditions. 


At  resonance 


(55) 
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Eq.  36  gives  the  mean  square  displacement  response  of  the  single  degree 
of  freedom  system  to  a  force  spectral  density.  We  shall  use  the  subscript, 
r,  to  denote  random  instead  of  the  subscript,  s,  which  previously  stood 
for  source.  Eq,  36  then  is 


If  we  define  the  applied  power  spectrum  per  1  cps  band  width  instead  of 
per  radian  per  second  we  have 


W 

r 


and 


Or 


S>=  wj! 


r  (*nj  | 


[221 

Fo 
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If  we  replace  yo/Fo  by  the  static  stress  per  unit  pressure,  S0,  and  change 
the  applied  spectral  density  to  the  square  of  the  rms  sound  pressure  in  a 
one  cycleband  width  we  obtain  the  square  of  the  rms  stress  response. 


<s3>~  p  2  s„2  *K 


48 


The  ratio  of  random  to  sinusoidal  stress  is 


(56) 


2  2  0 
S  =  rtf  5  "  p~ 
r  n  s  r 

0  2  e  2 
S  &  p 

s  r  s 
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If  the  damping  in  the  two  environments  is  the  same 


2  2 
S  »  *f  5  p 
r  n  rr 


(58) 


Smith  and  Malme  (Ref.  53)  point  out  that  the  acoustic  damping  in  the 
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siren  test  duct  may  not  be  the  same  as  that  for  a  structure  in  the  open. 

In  fact  they  preferred  to  use  the  actual  measured  response  to  random 
pressure  rather  than  relying  on  a  calculation  such  as  Eq.  56  which  includes 
the  damping. 


5B,  Correction  for  Non-Linearity 

Aircraft  panels  exhibit  nonlinear  stress  response.  Conventional  skin 
stringer  construction  is  nonlinear  because  of  the  membrane  stresses  which 
become  important  at  relatively  low  loads.  Honeycomb  structures  on  the 
other  hand  have  almost  linear  response. 

The  correction  for  non-linearity  would  not  be  necessary  if  the  siren 
test  were  conducted  at  the  peak  damage  stress  S^.  The  actual  test  is  run 
at  several  pressure  levels  so  that  a  plot  of  stress  versus  pressure  level 
can  be  made.  The  correction  for  non-linearity  to  be  applied  to  Eq.  58  is 
simply  to  multiply  Eq.  58  by 


2 

(59) 

where  the  subscripts,  pd  and  t,  denote  peak  damage  and  test  values  at 
failure.  Ordinarily  the  actual  stresses  at  the  location  of  the  failure 
will  not  be  measured  but  the  test  should  give  a  curve  of  relative  stresses 
versus  pressure. 

=  ft_  (60) 

S 

t 


where  the  primes  indicate  relative  stresses.  \  then  can  be  written 


can 


The  teat  pressure,  p^,  is  known  so  that  the  relative  test  stress, 
be  obtained  from  the  load  curve.  S  ^  is  obtained  from  the  peak  damage 
curve  for  the  desired  life  and  Sc  is  obtained  from  the  original  S-N  curve 
at  the  test  life.  Then  Eq.  60  is  used  to  calculate  .  Finally  p  ^ 

can  be  read  from  the  load  curve  at  and  Eq.  59  used  to  calculate 


5C*  Correction  for  Multimode  Response 

Although  it  has  been  pointed  out  that  only  a  few  of  the  many  possible 
modes  of  a  structure  will  be  important  there  is  the  problem  of  deciding 
how  many  and  which  of  the  inodes  are  significant.  Ideally  there  should  be 
enough  strain  gages  properly  located  tc  indicate  the  maximum  stress  in 
each  mode.  The  modes  producing  the  highest  stresses  are  the  most  import¬ 
ant. 

Belcher  (Ref.  40)  stated  that  probably  the  greatest  errors  in  the 
sine  -  random  equivalence  computation  were  encountered  in  the  interpretation 
of  multiple-mode  data.  In  multi-mode  response  the  stress  response  cycles 
are  not  fully  reversed  as  they  are  if  a  single  mode  response  is  predominant. 
The  suggested  multimode  correction  factor  to  be  applied  to  Eq.  58  is  the 
ratio  of  the  total  rms  stress  to  the  rms  stress  in  a  particular  mode.  A 
separate  correction  is  required  for  each  mode. 


The  basic  information  required  is  the  stress  versus  frequency  curve 

with  constant  pressure,  p  f  applied.  The  damping  ratio,  g,  and  the 

s 

relative  stress  for  each  mode  can  be  obtained  from  it.  Substituting  for 
2 

from  Eq.  58  we  obtain  for  constant  p  , 
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5D.  Corrected  Sine  -  Random  Stress  Equivalence 

The  corrections  for  nonlinearity  and  multi-mode  response  as  given  by 
Eqs.  61  and  62  applied  to  Eq.  58  yield 

Sr2  2 

2  "  *fn  6  £r  (63) 

9  V 

It  should  be  recalled  that  this  condition  applies  at  the  resonant  frequency, 

fn.  The  other  quantities  are: 

Sr  *  rms  8 tress  response  to  random  pressure 

Ss  *  rms  stress  response  to  sinusoidal  pressure 

p  =  random  rms  pressure  spectral  density 

(rms  pressure  in  a  one  cycle  per  second  band) 

p  B  rms  sinusoidal  pressure 
s 

8  ■  damping  ratio 

\  »  nonlinearity  correction,  Eq.  61 

Y  **  multimode  correction,  Eq.  62 

6.  Siren  Testing 

The  siren  testing  of  the  critical  structure  is  the  heart  of  the 
method  proposed  by  Belcher  (Ref.  40)  and  Fitch  (Ref.  39).  It  provides 
information  on  the  resonant  frequencies  and  shapes  of  the  important  modes. 
The  sample  structure  is  tested  to  failure.  The  life,  used  with  a  random 
S-N  curve,  in  effect  provides  the  correction  for  the  unknown  stress 

i 

concentrations  and  coupling  factors. 
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The  samples  to  be  tested  must  be  selected  so  that  they  are  represent¬ 
ative  of  the  critical  parts  of  the  aircraft*  This  depends  on  both  the 
distribution  of  the  sound  pressure  field  and  estimates  of  the  acoustic 
strength  of  the  various  parts  of  the  structure;  the  parts  with  the  highest 
ratios  of  pressure  to  strength  should  be  selected  for  study* 

It  is  of  course  the  objective  to  simulate  the  structure  and  its 
boundary  conditions  and  attachment  to  other  parts  of  the  aircraft*  This 
means  that  the  method  of  fastening  the  panel  should  be  the  same  as  in  the 
prototype.  Although  some  preliminary  work  may  be  confined  to  testing  of 
panels  only)  the  final  tests  should  include  other  parts  attached  to  the 
panels  which  might  influence  the  response  or  the  resonant  frequencies.  The 
electrical  hydraulic  and  other  hardware  must  be  included.  For  multiple 
span  panels  the  center  spans  are  more  likely  to  have  representative  boundary 
conditions.  Therefore  it  is  reasonable  to  strengthen  the  edge  spans 
slightly  by  decreasing  their  width.  This  will  increase  the  probability 
that  the  failure  will  occur  in  the  more  representative  spans. 

The  sample  should  be  adequately  strain  gaged.  Preliminary  frequency 

scans  using  a  low  pressure  source  outside  the  siren  test  chamber  are 

useful  for  determining  the  resonant  frequencies  and  observing  the  mode 

shapes*  Experience  with  such  tests  will  help  assure  that  strain  gages 

have  been  properly  placed  to  measure  the  maximum  stresses.  In  addition 

such  a  test  will  provide  the  data  on  damping  in  the  various  modes.  The 

frequencies  and  modes  that  produce  the  highest  relative  stresses  are 

chosen  for  the  siren  test.  If  more  than  one  mode  is  considered  important, 

tests  must  be  run  at  each  such  reasonant  frequency.  If  the  stress  versus 

frequency  curve  shows  nonlinearities,  such  as  unsymraetrical  resonance  peaks/ 

the  damping  must  be  determined  from  decay  curves  with  appropriate  filters 
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used  to  separate  the  modes. 

The  actual  fatigue  test  is  conducted  in  a  siren  test  facility.  Here 
grazing  incidence  is  used  because  it  usually  simulates  more  accurately 
the  service  environment.  Also  there  is  better  coupling  to  higher  order 
modes  of  the  structure  and  the  system  is  not  plagued  by  variations  of  the 
pressure  due  to  shifts  in  standing  ware  patterns. 

Stress  versus  pressure  data  are  needed  to  provide  corrections  for 
nonlinear  response.  This  can  be  obtained  by  the  technique  of  step  testing. 
Starting  at  some  relatively  low  sound  pressure  level  the  sample  is  exposed 
for  some  nominal  period  of  time,  say  15  minutes^  to  each  of  the  selected 
test  frequencies.  Then  the  sound  pressure  is  increased  a  small  step  and 
the  tests  at  each  frequency  are  repeated.  Fitch  (Ref.  39)  used  3  db  steps 
and  Smith  and  Malme  (Ref.  53)  used  2  db  steps.  The  entire  process  is  re¬ 
peated  until  failure  occurs.  The  test  results  give  the  cycles  applied 
at  each  frequency  at  each  pressure  level  as  well  as  maximum  relative  stress 
at  each  frequency  and  pressure  level.  This  technique  minimizes  the  number  of 
sample  which  must  be  tested. 

Fitch  (Ref.  39)  gives  the  following  discussion  summarizing  the  accuracy 
of  the  method  and  the  sources  of  error. 

"Comparison  of  test  resuits  under  random  and  sinusoidal  loading  has 
been  made  for  a  number  of  specimens.  The  variation  between  measured 
and  computed  stress  ratios  was  found  to  be  on  the  order  of  3  db.  Some 
of  the  more  obvious  sources  of  error  in  computations  for  stress,  and  for 
fatigue  life,  net  necessarily  the  order  of  importance,  are: 

1.  An  error  of  one  db  in  sound  pressure  measurement  represents 
approximately  12  per  cent  error  in  load. 
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2.  If  the  siren  excitation  frequency  is  off  resonance,  a  large 
non-conservative  error  in  damage  accumulation  can  occur* 

3.  Damping  factors  depend  on  how  they  are  measured* 

4.  The  propagation  direction  of  the  sound  relative  to  the  panel  in 

a  siren  test  and  in  an  air  frame  application  is  not,  in  general,  the  same. 

5.  Harmonics  of  the  siren  fundamental  pressure  wave  may  excite  higher 
modes  of  the  structure. 

6.  The  nonlinearity  of  the  structure  depends  not  only  on  the  design 
but  also  on  the  quality  of  fabrication,  which  is  variable  among  specimens, 
e.g.  skins  which  are  tightly  stretched  begin  to  diaphram  at  lower  pressures 
than  do  loose  skins.  This  can  have  a  large  effect  on 

7.  If  there  is  more  than  one  significant  mode,  additional  effects 
which  contribute  to  errors  exist. 

(a)  It  is  not  necessary  to  know  the  actual  values  of 
stress  for  each  mode,  but  the  relative  stress 
amplitudes  must  be  known  if  the  computed  value  of 
Y  is  to  be  meaningful. 

(b/  The  possibility  of  obtaining  misleading  strain  gage 
readings  because  of  a  non- zero  geometric  angle  between 
the  principal  stresses  must  be  considered. 

(c)  There  is  no  certainty  that  the  structural  area  which  is 
critical  when  all  modes  are  excited  simultaneously  (as 
by  random  noise)  is  the  location  of  failure  in  the 
discrete  frequency  test. 

(d)  Coupling  between  modes,  especially  when  there  is  little 

difference  between  the  resonance  frequencies,  causes 

difficulties  in  measuring  the  damping  factors  and  results 
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and  results  in  questionable  interpretation  of 
their  physical  meaning. 

8.  For  a  specified  life,  allowable  stress  varies  as  much  as  +  15 
per  cent  for  a  plain  smooth  specimen,  and  an  additional  variation  of 
+  15  per  cent  occurs  for  a  notched  specimen." 

7 .  Example  Test  and  Evaluation 

Fitch  (Ref.  39)  gives  an  analysis  of  a  hypothetical  problem  to 
illustrate  the  methods  he  had  proposed.  His  analysis  will  be  summarized 
here  for  the  same  reason. 

The  vehicle  is  assumed  to  be  a  Mach  3  intercept  fighter  having  two 
engines,  in  the  30,000-pound  thrust  class,  with  after  burners.  The  sound 
field  was  calculated  but  since  the  methods  were  illustrated  in  Section 
IV-2  for  the  A3J  aircraft  the  details  of  the  calculation  will  not  be 
repeated  here.  We  assume  that  we  have  the  overall  sound  pattern  shown  in 
Fig.  35.  The  specimen  to  be  simulated  and  fatigue  tested  was  selected  in 
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Detailed  analyses  of  the  varfrma  possible  -.lesions  of  the  aircraft 

were  made  to  estimate  the  relative  amounts  of  time  to  be  spent  at  the 

various  power  settings  of  the  engines  for  the  many  operations  performed. 

Detailed  operating  logs  for  various  types  of  aircraft  and  missions  are 

needed  to  accomplish  this.  The  key  information  derived  from  the  analyses 

for  our  present  purpose  was  that  there  would  be  an  estimated  22  hours  of 

operation  at  Max  A/B  power  at  sea  level  for  a  3000  hours  flight  life  of 
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figure  35.  Over-all  Near  Field  (Ref.  39) 


the  aircraft*  This  includes  take-^ff  and  maintainance  operations*  It  will 
be  recalled  that  the  engine  sound  pressure  level  decreases  rapidly  as  the 
aircraft  gains  speed  and  altitude. 

The  random  S-N  curve  for  the  6A1-4V  Ti  material  of  which  the  specimen 
was  constructed  is  shown  in  Fig*  36.  Th:  results  of  a  representative 
frequency  scan  at  143  db  is  shown  in  Fig*  37.  It  is  seen  that  two  resonant 
frequencies,  632  cps  and  695  cps,  appear  to  be  important.  The  specimen 
was  siren  tested  for  15  minute  intervals  at  each  of  the  two  resonant 
frequencies,  then  the  sound  pressure  level  was  raised  3  db  and  the  process 
repeated  until  failure  occurred.  The  stress  levels  were  recorded  at  all 
gages  during  the  runs*  The  test  results  for  the  gage  nearest  the  failure 
and  oriented  at  right  angles  to  the  line  of  failure  were  given  as  follows: 


Time 

Sound  Pressure 

Frequency 

Relative  Stress 

Minutes 

db 

cps 

psi 

15 

148 

652 

13,200 

15 

148 

695 

6,500 

15 

151 

652 

15,300 

15 

151 

695 

9,200 

3 

154 

652 

21,700 

The  plot  of  the  stress  load  curve  from  the  above  data  is  shown  in 
Fig.  38.  For  the  principal  mode  nt  652  cps  there  were  n  ^  *  5.9(10)"* 
cycles  at  151  db  and  =  1.2(10)"*  cycles  at  154  db.  The  condition  for 

failure  is 

<n151/N151)  +  (n154/N154)  =  1 


Assume  for  the  moment  that  all  the  damage  accumulated  at  154  db  so  that 
-  r‘154  =  1-2(10)“*.  From  the  S-N  curve.  Fig.  36,  the  corresponding 


stress  is  S 


154 


87,500  psi.  For  the  652  cps  mode  the  stress  load  curve, 


Fig.  38,  gives 


*151  =  -ill  X  Si54  =  1^122  *  S154  =  0-77  S 

e  t  oi  -Jin 


154 


21, 710 


154 


Then  S 


K1  =  67,300  psi.  The  S-N  curve  shovs  Chat  N, =  so  . 
151  -1U- 


This  means 


MATERIAL  6A1-4V  TITANIUM  Ffll  =  170,000  PSI  TYPE  OF  LOADING  REVERSED  BENDING 


(ISM)  SS3HJ.S 
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Figure  36  .  Calculated  Random  Fatigue  and  Peak  Damage  Curves  (Ref.  39) 
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Figure  37.  Stress  -  Frequency  Curve  (For  Example  Problem)  (Ref.  39) 
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RELATIVE  STRESS 


Figure  38.  Stress  -  Load  Curves  (For  Example  Problem)  (Ref.  39 
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that  there  was  no  damage  accumulated  at  151  db.  If  had  not  been 

infinite  it  would  have  beeu  necessary  to  solve  for  and  by  trial 

and  error. 

The  correction  for  multimode  response  is  Eq.  62. 
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The  frequency  scan  curve.  Fig.  37,  gives  the  following  data: 
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The  nonlinearity  correction  factor  from  Eq.  59  is: 

-.2 

Spd  ft 
St  PPd| 
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The  desired  life  of  the  structure  was  found  from  the  mission  analysis  to  be 


22  hours  or  5.2(10)  cycles  at  652  cps  for  Max  A/B  operation  at  sea  level. 
At  5.2(10)^  cycles  the  peak  damage  curve  in  Fig.  36  gives  Sp^  =  90,000  psi. 
At  the  test  life,  -  1.2(10)"*  cycles,  the  original  S-N  curve  in  Fig.  36 

gives  =  87,500  psi. 


From  Eq .  60  we  find: 

,  / 

*pd 


S/j  "  if  ^  *  IftH  x  21,710  -  22> 330  psl 


At  Sf,  =  22,330  psi  the  stress  load  curve  gives  p  =  0.151  psi  for  the 


pd 


pd 


652  cps  mode.  The  test  sound  pressure  level  is  equal  to  0.144  psi. 

.20,000  „  CKU4]2  =  0^963 


\  = 


87,500 


0.151 
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From  the  siren  random  equivalence  relation,  Eq.  63,  we  obtain 


jc&f 


n  / 


I  i 

1  rx  I 


f 


Here  »  87,500  psi  or  61,800  psi  (rms).  The  random  stress  “ 

23,000  psi  (rms)  is  found  from  the  random  S-N  curve,  Fig.  36,  at  the  desired 
life  of  5.2(10)^  cycies<  Then  with  8  “  0.0046,  f  ■  652  cps,  1.12, 

\  “  0,963  we  obtain 


—  =  0.117 

Ps 

db  =  db  +  20  log  _  0.117 

r  s  10 

db  =  154  -  19  =  135 

r 

This  is  the  allowable  random  spectrum  level  at  652  cps  for  the  desired  life. 
Vith  a  flat  spectrum  the  600/1200  cps  allowable  octave  band  level  would 
be 

^600/1200  *  +  *°®10  ^ctrave  width  =>  600  cps) 

=  135  +  28  »  163  db. 

The  acoustic  environment  at  the  location  simulated  by  the  test  specimen  was 
given  previously  with  150  db  in  the  600/1200  cps  band.  Thus  there  is  a 
13  db  margin  which  should  be  adequate. 

A  similar  calculation  for  the  695  cps  mode  indicated  a  design  margin 
of  8  db  with  the  assumption  that  failure  had  occured  in  this  mode  at  the 
end  of  the  151  db  test  period.  This  did  not  happen  so  we  know  the  design 
margin  must  be  greater  than  8  db.  In  order  to  find  out  more  accurately 
what  it  is  a  step  test  must  be  run  on  another  specimen  in  this  mode  only. 
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8.  Design  Charts 


Belcher  (Ref.  40)  and  more  recently  McGowan  (Ref.  41)  have  described 
the  development  of  design  charts  or  nomograms  which  relate  fatigue  life, 
sound  pressure  level  and  the  dimensional  parameters  of  several  of  the 
simpler  structures. 

With  a  J-71  jet  engine  as  the  source  of  random  noise  Belcher  and  his 
associates  obtained  a  limited  number  of  fatigue  failures  of  skin-rib 
constructions.  They  used  the  data  as  the  basis  in  developing  a  design  chart. 

For  a  long  panel  the  stress  in  the  skin  is  proportional  to  the  pressure 
and  the  square  of  the  ratio  of  width  to  thickness  (rib  spacing  to  thickness) 


8 


Where  o  =  rms  stress,  p  =  rms  pressure,  S  -  rib 

Using  this  for  the  ratio  of  stress  to  pressure 
that 


spacing, 


t 

s 


t  -  skin  thickness, 
s 

in  Eq.  56  we  see 
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Assuming  that  the  ratio  n/5  is  constant  for  panels  of  various  widths  and 
thicknesses  it  is  possible  to  draw  a  nomogram  relating  the  stress,  pressure, 
rib  spacing  and  rib  thickness.  Belcher  introduced  the  fatigue  life  data  from 
the  above  tests  in  the  same  design  chart.  Although  he  did  not  show  that 
^n/5  is  constant  the  data  appeared  to  fit  within  the  normal  scatter  for 
S-N  data.  Siren  fatigue  data  were  converted  by  the  sine- random  equivalence 
procedure  discussed  earlier  and  also  entered  on  the  design  chart  with  about 
the  same  accuracy  as  the  random  test  data. 

McGowan  evidently  carried  the  work  further  using  the  same  techniques. 

He  reported  that  the  ribs  were  critical  in  the  bend  radius  of  the  flange. 


Since  the  same  parameters  are  involved  for  the  ribs  as  for  the  skin,  the 
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rib  data  can  bt  superimposed  on  the  same  design  chart,  Fig.  39.  For  a 
long  panel  the  resonance  in  the  fundamental  mode  is  a  function  of  the 
8k4 n  thickness  and  rib  spacing.  This  too  is  shown  on  the  design  charts. 

The  design  chart  can  be  used  to  find  a  skin  thickness,  t^,  rib 

thickness,  t^,  and  rib  spacing,  S,  to  withstand  a  given  sound  pressure 

level,  db  ,  for  a  given  life,  N  *  Conversely,  given  the  panel  configuration 
K  R 

one  can  find  the  life  for  a  given  sound  pressure  level  or  the  allowable 
sound  pressure  level  for  a  desired  life*  Here  db^  is  the  rpectrum  level, 
db  per  one  cycle  per  second  band  width. 

Normally  the  life  is  given  in  hours  so  that  a  frequency  must  be 
estimated  to  obtain  the  corresponding  cycles  of  life.  The  sound  pressure 
level  for  the  frequency  is  obtained  from  a  plot  of  the  sound  level  spectrum 
and  a  safety  margin  added  to  give  db  *  The  design  chart  can  then  be  used 
to  obtain  several  combinations  of  rib  spacing  and  rib  and  skin  thicknesses 
to  meet  the  required  life.  After  a  practical  combination  has  been  selected, 
the  resonant  frequency  is  read  from  the  chart  and  a  new  cycles  of  life 
computed.  A  new  db  is  then  read  from  the  chart  and  compared  with  the 
sound  pressure  level  at  the  resonant  frequency*  If  the  design  margin  is 
inadequate,  some  other  combination  must  be  tried, 

McGowan  gave  design  charts  for  the  following  structures: 

a.  Conventional  skin  and  rib  construction. 

b.  Skin  and  rib  with  scalloped  doublers  at  the  ribs. 

c.  Attachments. 

d.  Rib  with  lightening  hole. 

e.  Beaded  panels. 

f.  Honeycomb  panels. 


RELATIVE  STRENGTH 


The  design  chart  for  the  skin  and  rib  with  scalloped  doublers  at  the 

ribs  is  given  in  Fig.  40.  The  doublers  improve  the  fatigue  characteristics 

9 

of  the  panels*  For  a  life  of  10  cycles  an  0.032  inch  skin  with  rib  spacings 
of  4  inches  will  withstand  123  db  without  doublers  but  will  withstand  131  db 
with  doublers.  The  skin  thickness  would  have  to  be  increased  from  0.032  to 
0.050  inches  to  accomplish  the  same  improvement  without  doublers.  Of 
course  the  rib  thickness  and  the  resonant  frequency  would  also  be  changed. 

The  design  charts  thus  facilitate  the  evaluation  of  the  design 
alternatives  for  panels.  This  is  very  helpful  because  most  of  the  failures 
are  panel  failures  near  or  aft  of  the  engine  exit  nozzle.  Although  the 
charts  have  been  adjusted  for  the  unknowns  such  as  stress  concentrations 
and  boundary  conditions  by  the  use  of  test  data  there  still  remain  ry  ,stions 
of  accuracy  of  simulation  and  the  errors  produced  by  multiple  mode  response. 
McGowan  concluded  that  the  charts  provide  the  means  to  meet  the  preliminary 
design  acoustic  requirements.  However,  if  the  charts  indicate  acoustic 
fatigue  is  the  primary  design  criterion  or  if  there  is  a  design  margin  less 
than  6  db,  McGowan  recommends  that  a  proof  test  be  conducted  in  a  simulated 
service  environment. 
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relative  strength 


ED 

FIGURE  40.  DESIGN  CHART -SKIN  WITH  BONDED  DOUBLER 
AND  RIB  CONSTRUCTION  (Ref.  41) 
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